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ABSTRACT 


This report, prepared by Martin Marietta Denver Aerospace, under Contract 
MAS9-I5302, presents the results of an analytical and experimental study of 
Space Shuttle propellant dynamics during External Tank/Orbiter separation in 
a return-to~launch-site (RTLS) abort situation. The study vas performed from 
April 1977 to January 1981 and was administered by the National Aeronautics 
and Space Administration, Lyndon B. Johnson Space Center, Houston, Texas, un- 
der the direction of Mr. Mark Craig, Technical Monitor. 
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1.0 


INTRODUCTION 


The Space Shuttle has been designated as America's prime space launch 
vehicle for the 1980s and beyond. The Shuttle system is a manned 
flight system requiring extensive mission planning, including intact 
operations. One such abort mode requires the mated orbiter/external 
tank (ET) to "fly” back to the launch site st an altitude of over 61 km 
(200,000 ft) using the main orbiter engines. This return to launch 
site (RTLS) mission abort sequence requires the orbiter and ET to sepa- 
rate under aerodynamic loading when a significant amount of propellant 
remains in the ET. The typical separation sequence is as follows: 

1) Begin powered pitchdown, from angle-of-attack >40 deg, at MECO 
-20 s; 

2) Achieve an angle-of-attack of -4 deg, using main engine thrust vec- 
tor control and the reaction control system; 

3) MECO, main orbiter engine cutoff; 

4) Coast for approximately 13 s using the reaction control system 
(RCS) to maintain attitude; 

5) Separation of the orbiter and ET using all downfiring orbiter RCS 
thrusters to move the orbiter away from the ET. 

During the separation sequence, the ET nominally contains a 1% residual 
volume of liquid oxygen (LOX)*; approximately 5.52 m3 (195 ft3), 

6,067 kg (13,376 lb). Additionally, the LOX line is full. The ET liq- 
uid hydrogen (LH 2 ) tank also contains a residual volume, but its im- 
pact due to liquid motion is small in comparison to the LOX tank, due 
to the low density of LH 2 . 

Aerodynamic forces decelerate the orbiter/ET during the 13-s coast 
after MECO. During this coast, the residual propellant "sloshes" for- 
ward developing a velocity relative to the ET at separation. Clearance 
between the ET and orbiter is provided, somewhat, by orbiter RCS jets 
moving the orbiter away, but more significantly, by aerodynamics acting 
on the ET to move it away. During separation, the motion of the pro- 
pellant in the ET significantly influences the ET motion, which could 
cause the ET to recontact the orbiter. It is not possible to lower the 
dynamic pressure at separation sufficiently to avoid this phenomenon, 
nor can initial conditions be obtained to avoid it. 


*0ne percent is the current OFT 1 baseline volume. Future missions 
may baseline much higher separation volumes (e.g., 25X) . 
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NASA simulations of orbiter/ET RTLS separation indicate a potential 
collision problem. However, these simulations have been performed us- 
ing a propellant slosh model, which has not been verified by test. To 
validate these simulations, a need to obtain empirical data was identi- 
fied. The objectives of this study were: (1) to obtain test data in a 

low-g environment on scale model tanks for use in validating a slosh 
simulation approach; and (2) to develop a simulation approach and com- 
puter code for use by NASA in RTLS simulations. To achieve these ob- 
jef r've's, a three-phase study was initiated in April 1977. 

Ph'ijif! I of this study was a drop tower test program designed to simu- 
tin ET LOX motion during RTLS separation. Thirty-two tests were 
conducted in Martin Marietta's Drop Tower Test Facility using two 
1 /60t^h~sc£>le models of the ET LOX tank; one with internal baffles, 
and one without. A scaling approach developed in a previous NASA study 
[Kefs i, 2] was used to relate drop tower test time to full scale ET 
rime, and to insure liquid motion representative of the full-scale ET. 

During the Phase I tests, small biaxial accelerations were applied to 
the model tanks, simulating aerodynamic deceleration of the ET during 
the RTLS separation sequence. The resulting propellant reorientation 
forces exerted on the tank were measured by crystal load cells, and 
liquid motion was photographed at 200 frames/s. The tests were conduc- 
ted both with and without LOX tank baffles to facilitate analytical 
model development and to assess the effect of baffles on reorienta- 
tion. In addition, a limited number of tests were performed simulating 
inflow to the tank from the LOX feedline. Subsequent analyses led to 
the conclusion that inflow would be negligible during separation. 

Based on scaling considerations, three test liquids were chosen for use 
in the test program: FC114B2, FC43, and HEXANE. The Phase I study is 

detailed in Reference [3]. 

The Phase I study indicated the importance that Bond and Reynolds num- 
ber scaling had on the character of propellant reorientation. Tests 
demonstrated that the motion of the bulk liquid did not change, even 
though the Reynolds number was varied an order of magnitude by changing 
test liquids. However, both the Bond and Reynolds numbers for the mod- 
el tanks, were much smaller than the values for the full-scale ET LOX 
tank. Closer approximation of full scale Bond and Reynolds numbers was 
not possible due to restrictions in model size imposed by the drop 
tower geometry. 

To further investigate the effects of Bond and Reynolds number and to 
validate the scaling concepts used, it was desirable to conduct a test 
program with a larger LOX tank model. Hence, Phase II of the study was 
initiated. In this phase, an available l/lOtb-gcale LOX tank model 
was modified for use in NASA's KC-135 "zero-g" test aircraft. Aircraft 
testing allowed larger accelerations to be applied to the model tank. 
This resulted in a closer simulation of full scale Bond and Reynolds 
numbers, as indicated in Table 1-1. 


1-2 



ORIGINAL PAG” B 
OF POOR QUALITY 

TABLE 1-1. 

COMPARISON OF TEST AND FULL-SCALE PARAMETERS FOR PROPELLANT 
REORIENTATION TESTS 


Parameter 

Full-Scale 

LOX 

1/lOth Scale 
Freon 113 

l/60th Scale 
Freon 114B2 

Bond Number 
Reynolds Number 

2.19 X 105 
5.16 X 107 

2.70 X 10^ 
2.39 X 106 

4.49 X 102 
1.35 X 105 


The Phase II test program consisted of four KC-135 flights involving 89 
low-g parabolas with an average test time of 20 s. Tests were conduct- 
ed both with and without the slosh baffles included in the l/lO^h- 
scale tank. Two test liquids were used to evaluate viscosity effects: 
FC113, and a water-methocel solution. The tank was mounted in an in- 
strumented fixture that measured the reorientation forces. In addi- 
tirn, liquid motion was photographed at 60 frames per second (fps) 
Adjustment of the aircraft flight path allowed application of biaxial 
accelerations to the tank, similar to those used in the Phase I drop 
tower study. Detailed documentation on the Phase II study is contained 
in Reference [4]. 

After evaluating the Phase I and Phase II test data, it was concluded 
that the scaling approach employed is valid, and that the liquid motion 
observed in the drop tower and KC-135 tests is representative of that 
expected in the full-scale ET [Ref 4]. The character of the liquid mo- 
tion observed in the two test programs was the same even though there 
was a factor of six difference in geometric scaling. The differences 
observed were primarily due to the reduced influence of surface tension 
forces (i.e., higher Bond number) in the KC-135 tests, which led to 
more breakup of the liquid surface. 

Test data from both test programs (Phasr- I and Phase IJ' were correlat- 
ed to analytical simulations of the test conditions, using a "slosh" 
model similar to NASA's Space Vehicle Dynamic Simulation (SVDS) model. 
The analytical model employed in these correlation studies, LAMPS [Refs 
1, 2], portrays the propellant as a point mass moving on an ellipsoidal 
constraint surface. The constraint surface is a distorted ellipse that 
best fits the center-of-mass envelope obtain??, by analytically rotating 
the tank in a 1-g field. Results of this correlation study indicated 
that the point mass model is inadequate to represent the propellant in- 
teraction forces during reorientation. The model was found to be con- 
servative in predicting the reorientation forces. Both the magnitude 
and impulse (integrated effect) of the liquid forces predicted by the 
analysis exceed those measured in test. 

It was therefore concluded that the NASA SVDS slosh model will over 
predict the slosh forces acting on the ET during RTLS sen;»ration. 

Phase III of the study was initiated based on these resv s. 
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Phase III consisted of the development of two and three dimensional, 
multiple mass models for simulating Idw g propellant slosh. The model- 
ing concept developed, employs a finite element representation of the 
liquid. This formulation results in a realistic distribution of liquid 
mass within the tank. Finite elements are used to express liquid in- 
compressibility by maintaining a constant elemental volume. The ele- 
ments are also used in the formulation of internal viscous shear for- 
ces. The tank is represented by a ellipsoidal surface, and liquid 
"nodes" are constrained from penetrating this boundary. The formula- 
tion of this modeling approach and the assoc, ted computer codes are 
documented in this report. The models developed in Phase III have been 
used to simulate Phase I and Phase II test conditions and the results 
correlate very well to the measured data. 

This report summarizes the entire 3-year study. Chapter 2.0 presents a 
discussion of both test programs. More detailed discussion of the test 
programs can be found in References [3 and 4]. Chapter 3.0 details the 
two dimensional finite element model formulation [LAMPS2] and correla- 
tion. Chapter 4.0 details the three dimensional model formulation 
[LAMPS3] and correlation. Conclus'ons of the study and recommendations 
for further development are presented in Chapter 5.0. 
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SCALE MODEL TESTING 
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The primary objective of the scale model testing was to acquire em- 
pirical data, on the character of propellant reorientation during the 
RTLS orbiter/ET separation sequence, for use in validating an 
analytical simulation approach. The full-scale conditions were 
scaled so that representative liquid (LOX) motion could be produced 
in subscale tanks. Phase I testing used a l/60th scale ET LOX tank 
and Martin Marietta's Drop Tower Test Facility. Phase II testing 
used a 1/lOth-scale ET LOX tank and NASA's KC-135 "zero g" test air- 
craft [Refs 5, 61. 

i.'uring both test programs, liquid reorientation was photographed at 
high speed and the dynamic interaction forces applied to the tank by 
the moving propellant were measured. Comparison of data from the two 
test programs permitted validation of the scaling approach. This 
chapter summarizes both phases of the test program. More detail on 
the experimental aspects of this study can be found in References [3 
and 4] . 


2.1 SCALING ANALYSIS 


During the RTLS separation sequence, the ex irnal tank (ET) experi- 
ences axial I'x | and lateral (Z ) accelerat. >ns due to aerodynam- 
ic drag. ThesS accelerations a?e shown in Figure 2-1, with respect 
to the orbiter coordinate system. The initial position of the re>^i 
ual LOX is established by the direction of the main engine thrust 
vector just before MECO. The initial position of the LOX free sur- 
face is perpendicular to the thrust vector, which is oriented at an 
angle y to the orbiter X axis. The following values for each of 
these variables defined the iull-scale conditions that were consider- 
ed in the test program: 

x‘o - 0.015 to 0.030 g 

’z'o - -0.005 to -0.030 4 

y =• Oo, 13U, 30 ^ 

Residual LOX volume ^ l%„ IX, 5X, lOX, 15%, 25% 
at separation. 

A dimensional analy-is of variables related to liquid reorientation 
in a container yields ttie following dimensionless parameters that 
characterize motion: 

„ V (Frouric number); ratio of inertia , 

r r = — .. j _ - ( 2-1 

to gravity corces 
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go » (Bond number) ; ratio of gravity 

a to surface tension forces 


pVr (Reynolds number); ratio of Inertia 
U to viscous forces. 


(2-3) 


where 


A * Liquid acceleration (cm/s2) 

V “ Liquid velocity (cm/s) 

r * Characteristic dimension, tank radius (cm) 
p = Liquid density (gm/cm^) 
cf = Liquid surface tension (dynes/cm) 
p = Liquid viscosity (cp) 

With respect to large amplitude liquid motion, it has been shown iRef 
7] that the Froude number can be related to the Bond and Reynolds 
numbers. Based on numerous liquid reorientation tests, an empirical 
expression has been established to describe this relationship [Ref 7]: 


( r. 

Bo ) 


1 0.48 [] 

^_(0 84j4.7] 

(2-4) 

' Bo ' J ) 



Experimental data showing the relationship between Froude number and 
Bond number is shown in Figure 2-2a. These data indicate that the 
Froude number is independent of Bond number when the Bond Number is> 
10. However, adjustments for viscosity must be made at low Reynolds 
numbers (Fig. 2-2b). If the Reynolds number is greater than 50, vis- 
cous effects are negligible. Therefore, these relationships indicate 
that for any propellant reorientation that has a Bond number greater 
than 10 and a Reynolds number greater than 30, scaling can be based 
on Froude number alone. As will be shown later, the requirements for 
Froude number scaling are satisfied for both the full-scale RTLS sep- 
aration event and the scale model tests. 

Equating Froude numbers for the full-scale ET LOX tank and the scale 
model tanks allows calculation of the full-scale time simulated by 
the test. 


Fr = Fr 
p m 


(2-5) 
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The subscript "p" refers to the full-scale tank, and the subscript 
"tn" refers to the model. Therefore, equating Froude numbers: 

V V 

P_ = m (2-6) 

Ja r Ia r 

‘N p p “N m m 

Since velocity is proportional to the product of acceleration and 
t ime (V » At ) : 


t 

P 


t 

m 




A r 

A r 
p m 


, seconds 


(2-7) 


Equation (2-7) defines the time scaling for a selected geometrical 
scaling and the ratio of prototype (full scale) to model accelera- 
tions. It is independent of the liquid properties. 


2.1.1 Geometric Scaling 

Selection of model scaling was based on being able to simulate enough 
full-scale time to encompass the critical period during the RTLS sep- 
aration sequence. This critical period was defined by NASA to be the 
time during which the propellant moves from its initial position to 
the front of the ET LOX tank. The required simulation time is ap- 
proximately equal to the time required for an object to free fall the 
tank length (~15 m). Based on the expected RTLS accelerations, the 
desired full scale simulation time varies from - 10 s (for 0.03 g) to 
~14 s (for 0.015 g). 

The scaling analyses for both the Phase I (drop tower) and Phase II 
(KC-135) test programs had to consider facility-imposed limits on the 
values of r^, A^,, and t„,. A primary consideration was to ma;<i- 
mize the reorientation forces (to facilitate measurement) v»hile main- 
taining adequate full scale simulation time, tp. To maximize 
forces, a large model size was desirable. For both test programs, 
the maximum model size permitted by the facility was selected as op- 
timnm. For drop tower, a l/60th-scale ( r^/rp) monel was chosen, 

and for the KC-i.35 a 1/lOth-scale model was chosen. An additional 
consideration was that a 1/lOth-scale model was already in existence 
that could be modified for the test program. 

To maximize the time scaling and reorientation forces, a large value 
Ajn is desirable; however, some practical limits exist. In the 
drop tower, large values of A^ reduce the test time, tm, because 
the drop capsule is accelerated with respect to a drag shield and the 
travel distance is fixed. A value of A^, =» 0.09 g is a practical 
upper limit that will still provide adequate test time (-1.6 s) . 

In the KC-135, Aj^^ is limited to = 0.2 g by operational considera- 
tions (Ajjj is actually a negative acceleration resulting in adverse 
wing loading on the KC-135 |. 
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Based on model sizes chosen and test acceleration limits, the full- 

scale time can be calculated. Table 2-1 tabulates the full-scale 1 

time simulated for each second of test time in the drop tower and | 

KC-135 test programs. The tabulated values indicate that scaling is 
adequate to simulate the critical RTLS period | tp ~ 10 to 14 s| . 

Table 2-1. TIME SCALING FOR THE DROP TOWER AND KC-1S5 TEST PROGRAMS '' 


Test Program 


V ^8^ 

t /t * 
p' ‘'m 

Ap - 0.015 g 

Ap » 0.030 g 

Drop Tower 

60 

0.09 

18.97 

13.42 

KC-135 

10 

0.1 

8.16 

5.77 

KC-135 

10 

0.2 

11.55 

8.16 


*Approxiraate values of t^ are as follows: 
Drop Tower: 1.6 s 

KC-135: t^ ~ 20.0 s 


2.1.2 Test Liquid Selection 

The properties of the test liquids enter into the scaling by assur''‘g 
that Bond and Reynolds numbers are sufficiently large to meet the re- 
quirements for Froude number scaling. A very dense liquid helps make 
both the Reynolds and Bond numbers large, and also assures that the 
forces due to a given volume of liquid will be maximized. Low sur- 
face tension and viscosity are also desirable. Safety and compati- 
bility with the plastic tank models were other considerations in 
choosing the test liquids. 

For both test programs, several test liquids were chosen to evaluate 
the effects of viscosity. In the Phase I drop tower testing, three 
liquids were used: FC*-43, FC-114B2, and HEXANE. Due to its ideal 

properties, the primary test liquid for this phase was FC-114B2. The 
other two liquids were used in a limited number of tests to evaluate 
viscous effects. 

Freon 113 (also known as Freon TF) was the primary liquid used in the 
Phase n KC-135 testing. Freon 113 has a high density, low surface 
tension and viscosity, and does not pose any hazards. Freon 113 has 
a vapor pressure of 3.7 N/cm2 (5.3 psia) at 2QOC, which is 
usually not a problem at normal temperatures. However, problems of 
rapid vaporization were encountered when the Freon and model tank 
were at temperatures around 30®C. 


*FC = Fluorocarbon 


- Jk 
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To •valuato tha Influaaca of viaeoua af facta* a aacond Fliaaa II teat 
liquid with a hi^ar viaeoaity waa daairad. Thia waa achieved by 
uaing water to which a thickening agent bad bean added. Hethocel, a 
■ethyl calluloaa thickener aada by Dow Chanical* waa uaad to incraaae 
tha water viaeoaity. Type F4N Nathocal waa uaad in tha proportion of 
2.1 grawa par liter of water to incraaae tha viaeoaity from ita nor> 
nal 1.0 to S.l ep. Thia gave a factor of 10 incraaae in tha kinema- 
tic viaeoaity over tha Freon. An antifoam emulaion (SWS-211 made by 
SWS Silieonaa Corp.) waa added to the wixture to reduce foaming. 

The final criterion for teat liquid aelection waa that the Bond num- 
ber and Reynolda nuadier be aufficiently large to permit Proude num- 
bera acaling. Table 2-2 delinaatea the liquid parametera for the 
full-acale BT and the Fhaae I and Phaae II teat programa. 


TABLE 2-2. FULL SCALE AND TEST LIQUID PARAMETERS 


System 

Density 

(gm/cm') 

Surface 

Tension 

(dynes/cm) 

Viscosity 

(cp) 

Acceleration 

(r) 

Bo 

(^) 

Re 

LOX (-ISJ’C) 
Full-Scale El 

1.14 

1J.5 

0.195 

0.015 

0,030 

2.19 X 10^ 
4.39 X 10'’ 

5.17 X 
7.31 X 

10' 

lo” 

Phase I : 
l/60th Scale 
Drop Tower 
Tests 








- FC-ll'iB: 
I20°c;) ■ 0) 

.Mb 

18.8 

0.75 

0.09 

4.49 X )0‘ 

1.35 X 

10- 

- FC-Ai 

(20°t:) (11 1 

1.91 

In. 7 

b.50 

0 . 09 

4.91 X 10- 

1.3b X 

10** 

- HEXANE 

(jo'i;) 

0 . bb 

18.4 

0. a 

0.09 

1.54 X 10- 

9.30 X 

10" 

Phase II: 
1/lOth Scali 
KOIJ5 Tests 








- Fi;-ll I 

1.37 

JO. 1 

0.b7 

0.1 

0.2 

1.35 X 10" 
2.70 X 10" 

l.b9 X 
2.39 X 

10'’ 

10*' 

- »ater/ 
-Metiioce 1* 

0.998 

47.9 

5.05 

0.1 

O.J 

3.b0 X 10-^ 
7.20 X lO-^ 

1.42 X 

2.01 X 

10^ 

10-’ 

*Tj calculate Ke, the value to be used for velocity, V, must be defined. A 
representative velocity is the free-fall velocity based on tank length, which 
ai)pro.\ imates ihe maximum velocity achieved by the liquid as it first moves 

Che length of the tank. 


j ■^Pn>pe»"'. ie.s were me^isurod usioR actual samples 

of the test liquids. 
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As shown by th« valuta in Table 2'*2t tha basic raquirasMnCa for 
Froude number scaling were met for both teat programa. In addition* 
Che magnitudes of Bond and Reynolds numbers for FC-113 in the KC^ISS 
test program were within an order of magnitude of the values for tha 
fu'l'scale BT. 

The primary objective of the Phase II test program was to allow clos- 
er approximation of full-scale conditions, to validate the scaling 
approach. The dimensionless parameters for the aircraft testa t»ere 
significantly greater (closer to full scale) Chan those for the drop 
tower tests. Phase ll testing increased Bond number by approximately 
two orders of magnitude and Reynolds number by one order of magni- 
tude. As Bond number increases, the effects of surface tension de- 
crease, resulting in more surface breakup. 

In the drop Cower tests, Che Bond number was low enough to produce 
some surface cohesiveness, which was far less evident in the KC-135 
tests. When reorientation occurs at large Bond numbers (>103), the 
liquid motion may take the form of "rain" (primarily during axisym- 
metric reorientation). However, with off axis acceleration compo- 
nents, such as those present during RTLS separation, the lateral ac- 
celeration component acts to hold the liquid mass together as it 
traverses the tank wall. This type of notion has been demonstrated 
at the Bond number up to 6 x 103 under both axisymmeCric and off 
axis reorientation conditions [Ref 12]. 

In summary, scaling considerations and subsequent comparisons of drop 
tower and KC-133 test results, indicate that the test results should 
be representative of liquid reorientation during full scale RTLS sep- 
arat ion. 


TEST SYSTEM DESCRIPTION 


Phase I Teat Program 

The Phase I test program employed Martin Marietta's Drop Tower Test 
Facility. This facility was developed under numerous prior test pro- 
grams in order to simulate low-g environments. Only minor molifica- 
tions were required to perform l^aae I testing. 

Teat Hardware - The Martin Marietta drop tower test system consists 
of a drag shield that contains a drop capsule on which a teat module 
is mounted. Figure 2-3 is a schematic showing the major components 
of the drop test system. Figure 2-A is a photograph of the drop test 
facility. The system provides a maximum of~2 a of test time. The 
drag shield is released from a height ot 23 m (=75 ft), and is decel- 
erated by landing in a bin of wheat at the end of the free fall 
period. 
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Figure 2-5 is a photograph of the test module that was used in the 
Phase I test program. The test module is shown mounted on the drop 
capsule. It consisted of a model tank, force measurement module, and 
a slider mechanism. Figure 2-6 is a closer view of the measurement 
module in which the model tanks were supported. The measurement mod- 
ule supported the tanks on three crystal load cells, and allowed for 
tank rotation to achieve various initial liquid orientations. The 
calibration fixture (Fig. 2-6) was used to calibrate the load cell 
system. The load cells (Kistler model 912) were fully capable of 
measuring small reorientation forces (on the order of 0.4 N, 0.1 lb) 
and surviving impact decelerations at the end of the test. 

Two l/60th-scale replicas of the full scale ET LOX tank were con- 
structed; one with no internal baffles (Fig. 2-7), the other with 
both slosh and antivortex baffles (Figs. 2-8 and 2-9). Although the 
baffles were simplified somewhat, compared to the full-scale LOX 
tank, the overall geometry was accurate. Spacing between the baffles 
and tank wall was generally representative of the actual LOX tank. 
However, •’he model did not reflect the "zero gap" between the wall 
and baffles that occurs at two of the baffle stations in the full- 
scale tank. This had some effect on reorientation for very small 
fill volumes, allowing the liquid to flow primarily behind the 
baffles . 

Biaxial accelerations were applied to the model tanks by using linear 
spring motors that: (1) accelerated the measurement module laterally 

on the slider, and (2) accelerated the drop capsule axially with re- 
spect to the drag shield (Fig. 2-3). The drop capsule was automati- 
cally released from the drag shield at drop initiation. Three force 
links (containing the load cells) allowed forces within the plane of 
the applied accelerations to be measured. Bearings at each end of 
the links permitted only forces along the link axis to be measured. 
Three flexures, perpendicular to the plane of the force links, pre- 
vented any out-of-plane tank motion. These flexures were sized (by 
stiffness) such that they would not significantly alter the forces 
measured by the load cells. The design of the test system resulted 
in essentially planar reorientation. Out-of-planc motion and forces 
were insignificant. 

As mentioned in Chapter 1.0, some tests included the simulation of 
inflow into the tank from the LOX line. This was achieved by using a 
spring-operated piston/reservoir controlled by a solenoid valv . 

Flow rate was controlled by a metering valve. The tests inflowed a 
volume of liquid approximately equivalent to totally draining the LOX 
feedline. However, an analysis (Ref 3] has indicated that, conserva- 
tively, only about 5% of the liquid in the line will drain into the 
tank during RTLS. This volume (»0.05% of the LOX tank volume) should 
have a negligible effect on propellant reorientation. 

The liquid motion was filmed using a 16 ram Milliken DBM-3a camera, 
mounted on the slider, at a frame rate of 200 frames per second 
(fps). Test accelerations were varied by changing the magnitude of 
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the spring motors. The lateral acceleration applied to the tank was 
calculated based on lateral travel as a function of time. Lateral 
travel was determined by a scale mounted in view of the camera and 
time was based on camera frame rate. Average axial acceleration was 
calculated by kinematic equations for the two bodies (drag shield, 
drop capsule); given body masses, applied forces, and travel distance 
of the drop capsule relative to the drag shield. The travel time of 
the drop capsule relative to the drag shield was based on the load 
cell time histories. 


Test Conditions - Test conditions were chosen to provide reorienta- 
tion representative of that expected during RTLS separation. Table 
2-3 presents the expected ET accelerations due to aerodynamics during 
flLS: (1) during mated coast, and (2) after orbiter/ET separation. 

2'able 2-3. EXPECTED ET RTLS ACCELERATIONS 



ET Acceleration 
(Orbiter Coordinates*) 

RTLS 



Phase 

(g) 

Zc (g) 

Mated Coast 

0.0*5 

-0.005 

After Separation 

0.03 

-0.03 

♦See Fi 'ure 2-1 


In the drop tower, lateral and axial accelerations were applied to 
the tanks. Figure 2-10 shows the relationship between orbiter coor- 
dinates and the test-applied accelerations. These accelerations are 
related to orbiter coordinates as follows. 


* Xq cos y + ZiQ sin y (2—8) 

^lat * ~2o cos y+ "Xg sinr (2-9) 

The desired ratio of Ag^ial to Aiat oan be calculated for the 
RTLS accelerations of Table 2-3 ( y assumed 13 deg) using equations 
(2-8) and (2-9). 


Mated coast: 


^xlal^^lat 


1.64 


After Separation: 


A . ,/A, . = 
axial lat 


0.62 


Due to drop tower limitations and the desire to maximize test times 
and force levels, use of these actual ratios was not feasible. A 
value of A^xial * 0.09 g is been found optimal from test time and 
force level consideratici.a in previous test programs. Choice of a 
value for Aiat was based on the desire for a test time (tm) of 
about 1.6 s (Table 2-1). Because of the available lateral slider 
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trav«l, UcaraL accalcrationa graatar than S 0.03 g will raduca tha 
taat tina balow 1.6 a. Conaaquant ly , tha najority of tha taata uaod 
tha ratio: 

\xial/^lat " 0.03 * 

Although this approach did not duplicata actual axpactad ratioa of 
acralaration componanta, it did produca liquid motion rapraaantativa 
of full-acala conditiona, which waa tha primary objactiva. Aa long 
aa thara ia a lataral componant of accalaration, tha liquid will mova 
along tha tank walla rathar than through tha tank cantar. Changing 
tha ratio of accalaraciona primarily changaa tha final aquilibrium 
poaition for tha liquid and tha rata of motion. Tha baaic charactar 
of tha motion will ramain tha aama. Tha taat raaulta aupportad thia 
approach. 

Tha taat matrix waa atructurad auch that tha affact of varioua param- 
atara could be avaluatad. Table 2*4 praaanta tha actual taat matrix, 
including calculated accalaration valuaa. 

Data Reduction > To tacilitata correlation with analytical raaulta, 
c# Tt <iata wara racordad and latar tranaformad to a tank cantarad 
coordinate ayatam (Fig. 2-10). Digit iaad data wara amoothad to ra- 
roova taat fixture and camara-inducad noiaa with a moving average dig- 
ital low paaa filter aat at 10 Ha. Ait invaat igat ion waa conducted to 
inaura that filtering did not reduce peak force lave la. 

Data interprat.at ion and correlation had to conaidar tha fact that 
forcea due to the tank and support atructura maaa wara imbedded in 
the load cell data along with liquid reorientation forcea. Another 
cona idarat ion wa.a tha large tranaient encountered at drop initiation 
whan tha tank tranaitioned from normal gravity to low g teat accaler- 
ationa. To overcome theae problema during analytical correlation, 
tha analytical aimulationa account for actual teat conditiona (in- 
cluding loada due to atructura I maaa). Reference (3) dataila the 
actual iivita reduction procaaa. 

2 Miaae II Taat Program 

NASA’a KC-13S "aero-g*' taat aircraft waa uaad in tha Phaae II teat 
program to produce low-g environmenta repraaantat iva of RTLS aepara- 
tion. The teat fixture and aaaociated hardware ware fabricated and 
aaaembled by Martin Marietta and NASA JSO personnel. Both organiaa- 
tiona participated in teat conduct and data reduction. The objective 
of the Phase 11 taat program waa to more closely simulate the Bond 
and Reynolds numbers expected during full-scale RTLS separation to 
validate the scaling approach of Section 2.1. 

Teat Hardware - The teal fixture waa a scaled up version of the Phase 
I drop tower maaauraroent nxidule. Tha fixture consisted of a box 
frame that waa suppv^rted from a base structure by three strain gaga 
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TABLE 2-4. DROP TOWER TEST MATRIX 










Parameter 


Test 

Acceleration (g) 


l.nx 

»at'l led 

i Kill 

of 


l.i ‘III id 

■\ixlal 

•'l.tt 

ir 

Out t low 

lank 

Interest 

1 

1 

0.090 

0.029 

13 

Yes 

No 

2 


2 

2 

0.093 

0.031 

13 

Yes 

No 

2 

LOX 

3 

1 

0.093 

0.031 

13 

Yes 

Yes 

2 

Outtlow 

4 

2 

0.092 

0.032 

13 

Ves 

Yes 



5 

HKH 


0.030 

13 

No 

KSK 

K9 


6 



0.029 

13 

No 


H9 


7 




13 

No 


B9 

Liquid 

8 



0.032 

13 

No 

■H 

■1 

Viscosity 

■■ 


^KSKeSH 

0.031 

13 

No 

No 

■9 



2 


0.032 

13 

No 

No 



WM 



0.034 

13 

No 

Yes 

^B 



2 

0.092 

0.033 

13 

No 

Yes 


BHHBH 


1 

0.088 


msm 

No 

bbki 

2 

HpBi^B 


2 

0.088 


Bil 

No 


L_i 


15 


IKRE9i 


KB 


HOB 

wmM 


16 




m 




Cliaiii’e 

17 


0.090 


13 



B9 

in Aiat 

18 


0.091 


13 

bbh 

9RK 

BE9 



IHHHI 

91 

mai&jmm 

BK 


Yes 




1 



Kfl 


\ es 


■ nil 

tfl 

1 

0.091 


13 


No 



Ksl 

1 

0.093 

0.032 

13 

■9B 

No 

15 


■w 

1 — 

HSfli 

Mint ilflIB 




2 



1 

1 ■ 





2 

>-30° 

19 

1 

Bin ;B 

■ilfl 

^^B 

HiH 


10 


26 

1 

0.088 

0.031 

wEm 


■DK 

10 


WSM 

1 





Yes 

2 



1 



0 


No 

2 

7-0° 

29 

1 

0.050 

0.038 

0 

No 

Yes 

10 

'^axial 

30 

1 

0.042 

0.032 

0 

No 

No 

10 


31 

1 

0.090 

0.032 

13 

Ves 

Yes 

2 


32 

3 

0.090 

0.032 

13 

No 

No 

10 

Special 









lest s 


# Liquid 1 ■= FC114B2 

2 - FC43 

3 = HEXANE 


* Lateral slider delayed u.i s 

A Tank outlet increased by factor of 2, flow rate reduced 
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Figure 2-10. Orientation of Tank Coordinate System 
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lo«d cells (Interface, Inc.). Tlte tank was mounted to the box frame 
at two pivots, allowing the tank to be tilted at 0, 13, and 30o (>), 

The load cells mea- ured forces in the plane of the tank parallel with 
the aircraft longitudinal axis. Rod end bearings were used to con- 
nect the load cells to the base and box frame. Four flexures pre- 
vented out-of-plane movement of the tank and box frame. These were 
sized, by stiffness, such that they would not significantly affect 
load cell measurements .The tank was a 1/lOth scale model of the ET 
LOX tank. It's overall length was 151 cm (4.95 ft) and the barrel 
section diameter was 84 ca (2.76 ft). The tank, which was 
fabricated of transparent acrylic plastic, was made in two sections 
joined with bolted flanges at the junction of the lower dome. 

The tank contained removable baffles. The antivortex baffle was made 
of transparent acrylic plastic. The slosh baffles were fabricated as 
an integral unit from aluminum. Both the slosh and rntivortex baf- 
fles were geometric scale replicas of the full-scale tank. However, 
the slosh baffles lacked details of the supports between baffles (un- 
important in the test program). Installation of the slosh baffles 
accurately represented the baf fle-to-wall spacing of the full-scale 
LOX tank. 


Figures 2-11 and 2-12 show the tank mounted in the fixture at the 13- 
and 30-deg orientations respectively. Figure 2-13 is a photograph of 
the test system including the instrumentation rack. 

During a flight, the volume of test liquid in the tank was varied nu- 
merous times. Transfer of the test liquid between the model tank and 
a supply tank was performed with an electrical pump. Figure 2-14 is 



Overboard 

Dump 


Figui>e 2-14. Plumbing Sahematio 
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a plumbing schematic of the teat system. Interchanging lines at the 
valves adjacent to the pump permitted the direction of flow to be re- 
versed for tank draining. A calibrated sight gage on the supply tank 
was used to determine the volume of liquid transferred to the tank. 

The plumbing system employed check valves [0.7 N/cm2 (I psi) crack- 
ing pressure] as breather and relief valves. These valves were in- 
tended to prevent the pressure differential between the tank and 
cabin from exceeding +0.7 N/cm2, as aircraft altitude varied. How- 
ever, on the second flight (with the baffles installed) when Freon 
113 (IX volume) was used as the test liquid, rapid vaporization oc- 
curred when the cooler Freon contacted the warm tank walls and baf- 
fles. Due to the relatively long test time (5:20 s), the Freon col- 
lected over the relief valve inlet (top of tank) and sufficient tank 
venting was not possible. This resulted in tank failure (Fig. 2-15). 

The tank was repaired and a large vent port was added in the side of 
the tank. A flexible rubber disk closed the port, but any increase 
in tank pressure caused it to open. This permitted essentially un- 
restricted release of vapor, since the liquid was not usually orient- 
ed near the port. This design functioned satisfactorily, preventing 
any further problems with tank pressure. 

As part of the repair, the flange on the upper part of the tank was 
replaced with an aluminum ring that was glued and bolted to the orig- 
inal plastic tank section. Figure 2-16 is a photograph of the re- 
paired tank. 

The motion of the liquid was recorded with 16 mm Milliken motion pic- 
ture cameras operating at 60 fps. Three camera positions were used 
for complete coverage. A three-axis accelerometer measured applied 
accelerations. Load cell data were output on a chart recorder for 
quicklook data evaluation and were also recorded on magnetic tape for 
subsequent data reduction. The acceleration time histories were also 
recorded. Synchronization between the film and recorded data was 
provided by a bilevel signal. When a push button was depressed, 
lamps were illuminated in the camera field of view, and the change in 
level of a dc electrical signal was tape recorded. Figure 2-17 de- 
tails the test configuration in the KC-135 . 

Test Conditions - Low-g accelerations were applied to the tank by 
slightly modifying the normal KC-135 parabolic flight path. During a 
typical test, the following flight profile was followed. The parabo- 
la was initiated by diving the KC-135 and then pulling up to achieve 
the desired airspeed and angle of attack for the parabola. During 
the pull up, the vertical acceleration, which was perpendicular to 
the aircraft floor (7,AC, Fig. 2-11), reached levels of +1.8 to +2.0 g. 




























ORlGfriAL PAGi IS 
OF POOR QUALITY 




m 


ORIGiMAL PAGE 
OF POOR QUALITY 



t'tjure L'-i/ 
2-23 


Figure 2-17. KC-1Z& Teat Configuratxcn 
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This acceleration acted to settle the liquid at the bottom of the 
tank, duplicating the propellant location in the BT before MS CO. As 
the aircraft entered the parabolic arc, there was a rapid transition 
(*2 s) from the high positive acceleration .. the negative lo^.r-g test 
acceleration. During the test program, accelerations of » -0.1 g and 
s-0.2 g '.^AO were used. The lateral component (XAC) was obtained by 
adjusting t.he throttles so that the aircraft continued to accelerate 
as it entered the parabola. 

During some tests, out-of-plane (YAC) accelerations occurred. Be- 
cause of the variability of the accelerations, approximately three 
teats were performed for each test condition and tests with out of 
plane (YAC) liquid motion or similar defects were discarded. 

Test conditions, similar to those in Phase 1, were selected to facil- 
itate comparison of results. Table 2-5 delineates the Phase II test 
matrix. The test program included 30 test conditions (a total of 89 
low-g parabolas). Four KC-135 flights were floim to accossnodate con- 
figuration changes: baffle installation and removal, and test liquid. 

Data Reduction - The data reduction process for Phase II was similar 
to that of Phase I. The analog flight tapes, including force and ac- 
celeration measurement3 , were digitized and converted to engineering 
units. Approximately 25 s of time encompassing each test were digit- 
ized. The load cell data were then resolved to the tank centered co- 
ordinate system (Fig. 2-10) for correlation with analyses. Refetence 
[4] details Che Phase II data reduction process. 


2.3 SUMMARY OF TEST RESULTS 


This section is a summary of Che test results from both Phase I and 
Phase II. Observations on the character of the liquid motion, as well 
as comparisons of measured force data are presented. More detail on 
the test results can be found in References [31 and [4]. 

In general, it was found Chat the drop Cower provides a more controlled 
test environment. Test conditions are repeatable and acceleration en- 
vironments can be accurately specified. Although liquid reorientation 
forces in the drop tower are small, they are measurable, and photo- 
graphic coverage is excellent. 
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Testing in the KC'-133| on the other hand, provides sone special prob- 
lems. Acceleration environments are quite variable and depend substan- 
tially on pilot technique. Because of the large transients at the be- 
ginning and end of a parabola, force resolution is still a problem. 
Instrumentation imist accommodate the large transients yet still be sen- 
sitive enough to measure the reorientation forces. The large transient 
accelerations (on the order of 2 g) dictate a "beefy" teat fixture that 
partially obscures photographic coverage. However, the KC-135 does 
provide for substantially longer full-scale simulation time than the 
drop tower. In addition. Bond and Reynolds numbers are much closer o 
those expected in the full-scale ET. 

A comparison of the KC-135 and drop tower test data (Ref 4l resulted in 
the conclusion that the general character of the liquid reorientation 
was the same in the two test programs. Although some minor differences 
were evident (e.g., more surface breakup in the KC-135 due to higher 
Bond number), the bulk motion of the liquid was similar. Reference [4l 
also presented a comparative analysis between the drop tower and KC-135 
data, which demonstrated the validity of the time scaling analysis of 
Section 2.1. 

Analytical correlations (Chap. 3.0) to both the KC-135 and drop tower 
test data, also validate the scaling approach. The analytical model is 
based on the premise that Froude number scaling is valid. 

The following paragraphs present some observations on the influence of 
selected parameters varied during the two test programs. Only drop 
tower force time histories are presented for comparison. The variabil- 
ity of KC-135 accelerations makes it difficult to directly compare 
force time histories. 


2.3.1 Unbaffled Tank 


Liquid reorientation in the unbaffled tank was well behaved. In all 
cases, the biaxial acceleration caused the liquid to reorient along 
the tank wall. Previous test programs [Ref 141 have demonstrated that 
the liquid will reorient along the tank wall, even for purely axial 
accelerations, if the tank is tilted as little as 1 deg with respect to 
the axial acceleration vector. During drop tower testing, a wave 
formed in the liquid center early in the teat, indicating a tendency 
for the liquid to traverse the tank center. This wave, generally 
called a Taylor instability, forms because the motion of the liquid 
along the tank wall lags the bulk liquid motion produced by accelera- 
tion forces. These instabilities, which were generally short lived, 
rapidly joined the flow along the wall. This instability was not as 
evident in the KC-135 tests due to the larger Bond number and the gen- 
erally larger lateral acceleration components. 

Figures 2-18 through 2-25 contain single frame photographs from the 
high speed drop tower and KC-135 movies. The discrete time of each 
frame, from test initiation, is noted. These figures will be discussed 
in the following paragraphs. 
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Effect of Fill Volume - Fill volume did not influence the general 
character of the' liquid motion. Figurea 2-18 through 2-20 contain ain- 
gle frame photographa from the high apeed drop tower s>oviea for fill 
volumea of 2, 15, and lOX, reapectively (FC-114B2). In general, a liq- 
uid film led the bulk liquid motion around the tank. In the caae of 
the 2Z fill volume, the majority of the liquid waa diatributed in the 
film. The leading edge of the liquid, for aome fill volumea, complete- 
ly circulated the tank. However, the laajority of the liquid began to 
decelerate and collect after pasaing through the tank 'a ogive dome. 

With a longer teat time, the liquid would have collected and eventually 
reached an equilibrium position as determined by the acceleration 
vector. 

The ogive dome of the ET LOX tank crapped the liquid as it passed 
through, especially for larger fill volumes. The general reorientation 
character was one of filling and then emptying Che dome. However, at 
the smaller 2Z volume, the liquid flowed smoothly through the dome. 

This reorientation character is different than Chat observed for tanks 
having hemispherical domes (Refs 1, 2]. Hemispherical domes allow the 
liquid to flow through smoothly. The entrapment tends to dissipate 
energy through turbulence. This observation is consistent with the 
results of other drop tower tests on cylindrical tank compartments per- 
formed at Martin Marietta. The sharp corners of these tanks signifi- 
cantly damped the liquid motion. 

Figures 2-26 and 2-27 are plots of the measured force time histories 
^^YT* f«nk coordinate system) for fill volumes ranging from 2 
to 15Z (identical test conditions). Larger fill volumes exhibit larger 
forces that peak sooner than smaller volumes. Smaller volumes must 
traverse a longer distance to reach the same relative point in the 
tank, which results in a peak force time delay. However, fill volume 
do**! not have an effect on the overall character of the measured forces. 

Effect of Liquid Viscosity - The scaling analysis indicated that 
the influence of viscosity during reorientation should be negligible if 
the Reynolds number is greater than 50. To investigate viscous ef- 
fects, the lOZ fill volume case was run with three different test liq- 
uids in the drop tower test program. Figures 2-20 through 2-22 show 
the reorientation time histories for FC-114B2 (U • 0.75 cp), FC-43 
(w* 6.50 cp), and HEXANE (u* 0.33 cp). Figure 2-23 shows a similar 
KC-135 test condition using FC-113 (U* 0.67 cp). The general charac- 
ter of the bulk liquid motion was the same in all four tests. The pri- 
mary observable difference was in the motion of the liquid film that 
leads the bulk liquid motion. A comparison of the 0.4-s frsmes in Fig- 
■ es 2-20 through 2-22, shows that the lower the viscosity, the higher 
film velocity. HEXANE leads, closely followed by the FC-114B2; the 
FC-43 film lags both. 

Viscosity also influenced how much liquid circulated completely around 
the tank during the test. However, the bulk motion was not signifi- 
cantly altered even though viscosity varied by a factor of 19.7 and the 
Reynolds number varied by a factor of 25.7 over the four tests. This 
supports the scaling premise that viscous forces have a negligible ef- 
fect on bulk liquid reorientation in the RTLS Reynolds number regime. 
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Pigur* 2-28 « tore* coatp«ri«on for drop towor t««(s 6, «nd 32 

(Piga. 2-20 CO 2-22). Absolut# msgnitud* cowpsrisons cannot ba aad* 
dua to Cha diffarant liquid danaitiaa. Howavar. any viscosity affacts 
would probably appear in tha phasing of force paaKa. PC-43 is over 
eight Cistas as viscous as PC-U4B2. and almost 20 tisMS as viscous as 
HEXANE . The results, howavar, are inconclusive. No clear trend is 
present in Pigure 2-28. Tha Pyt plot supplies the most information. 

Por Pi'-ll4B2, the negative Pyr peak occurs slightly ahead (in time) 
of the peak for PO-43. This may indicate chat tha less viscous 
PC-IUB2 is moving a little faster than the PC-43. In general, how- 
ever, viscosity does not seem to have an appreciable effect on the 
force time histories. 

Effect of Tan k Orientation - Several tests were run at tank orten- 
tsHons of 0 an<^ TO deg. In addition, the axial acceleration component 
was halved in soma tests. These changes did not significantly alter 
the overall liquid motion characteristics. Pigures 2-24 and 2-2' show 
drop tv^wer tests 30 lY» and 2' ( T» 30*^3, respectively. 

The KC-13' test program also included tests at Y* 0, 13, and 30 dag. 

No appreciable ditference in the liquid motion was observed in Chose 
tests either. 


2-JJ 




ORIGINAL PAQL* Y8 
OF POOR QI'AL1T\ 


Y - 13* 

Fill VoluiM - 2X 

T«»t Liquid ■ FC-11AB2 

0,031 

0.090 s 














ORIGINAL PAGE 
OF POOR QUALlT’i 


Fill VoluM • lOZ 

Test Liquid - FC-114B2 
0.030 g-^-i 


Figure 2~20. Drop Tower Test 5 , Unbaffled Tank 
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Figure 2-2S. KC-135 Teat 1.2.2, Unbaffled Tank 
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Figure 2-26. Drop Touer Teat 26, Unbaffled Tank 
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Figure 2-26. 

Uribaffled Tankt Drop Tower Liquid Reorientation Force Time 
Hietoriee: Fill Volume Comparieon (Teete 5 , 5 , and 21) 
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b) Zj, Foynse Time History 
Figure 2~27. 

IMbaffled Tankt Drop Tower Liquid Reorientation Forae lime 
Hietoriee: Fill Volume Comparison (Tests 5, 23, and 22) 
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h) Zy Force Time History 
Figure 2-28. 

Unbaj'fled Tavkt Drop Tower Liquid Reorientation Force Time 
BistoHee: Liquid Vieooeity Cempariem (Teete 6, «, and W 
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2.3.2 Baffled Tank 


Incorporation of baffles into the tanks resulted in extremely turbulent 
flow. Figures 2-29 through 2-34 contain single frame photos from the 
high speed KC-135 and drop tower movies for selected baffled tank 
tests. Although the general reorientation character was similar to the 
unbaffled tank, the baffles significantly altered the reorientation 
forces. 

For small fill volumes (2X), the antivortex baffle caused a portion of 
the liquid to traverse the tank interior and impact the wall near the 
ogive dome. This did not occur in the KC-135 tests because the XAC 
acceleration began •: act while the ZAC acceleration component was 

still positive. This tended to position the liquid away from the anti- 
vortex baffle. Construction of the l/60th-scale (drop tower) tank baf- 
fles allowed the liquid to flow between the baffles and the tank wall. 
The 1/lOth-scale (KC-135) tank baffles were more representative of the 
full scale ET LOX tank and prevented this wall flow. 

Figures 2-29 and 2-30 show reorientation for similar test conditions in 
the drop cower and KC-135. Although the higher Bond number in the 
KC-135 test resulted In more surface breakup, the general reorientation 
character was similar. The baffles in the KC-135 tank did a better job 
of containing the liquid in the upper part of the ogive dome (see Fig. 
2-30, T = 5.6 to 6.4 s). This was due to the smaller baffle/wall gap 
in the 1/lOth-scale tank and would be the expected behavior in the 
full-scale ET. 

Figures 2-31 and 2-32 show the influence of fill volume. The effect of 
the antivortex baffle on small fill volume (as discussed above) can be 
seen in Figure 2-31. Figures 2-33 and 2-34 show tests performed ct 
Y= 0 and 30°, respectively. As in the unbaffled tank tests, no ap- 
preciable difference in reorientation character was evident. 

Turbulence caused by the baffles generally reduced the liquid reori- 
entation forces. Figures 2-35 and 2-36 present ® comparison of drop 
tower reorientation force time histories for similar baffled and unbaf- 
fled tank test conditions. Similar results were noted for all test 
conditions. Figures 2-37 and 2-38 show drop tower test 16 and 18, in 
which the lateral (Ajj^t) acceleration component was delayed for 
=0.5 s after test initiation. Figure 2-39 presents a comparison of the 
force time histories for these two tests. Note that the baffles also 
reduced the force levels for this case. 
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a) Yy Fopoe Time History 



b) Zj, Force Time History 
Figure 2-35. 

Drop Tower Liquid Reorientation Foi'ce Time Histories: 
Baffled/ Imbaffled Comparison (Tests 6 and 8) 
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b) Zj, Force Time Hietory 
Figure 2-33. 

Drop Tower Liquid Reorientation Force Time Hietoriee: 
Baffled/ Unbaffled Comparieon (Teats 16 and 18) 
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3.1 INTRODUCTION 


The liquid motion during the RTLS maneuver has previously been analyzed 
by a single point mass mechanical analog developed for NASA's Space 
Vehicle Dynamics Simulation Program (SVDS). This analog portrays the 
liquid as a point mass moving on an ellipsoidal constraint surface 

+ ( 1 I ■ 1. This constraint surface is defined as the 

\ ao - aiz / ^ 

ellipsoid that best fits the locus of center of mass locations pre- 
scribed by analytically rotating the tank in a 1-g field. Therefore, 
the constraint surface is a function of tank geometry and fill volume. 
The forces that constrain the CM to the ellipsoidal surface re repre- 
sentative of the forces the liquid exerts on the ET tank. A similar 
model, called LAMPS, was developed at Martin Marietta. The LAMPS model 
portrays the liquid as a point mass moving on an ellipsoid constraint 
surface. Conceptually, the LAMPS model expresses the forces of the 
moving propellant as, 



F 



(3-1) 


where 


MA “ D'Alembert component of force; 


M 



modified centriptal component of force; 


f ~ effective mass factor, 

m 

The effective mass factor concept was developed in a previous NASA 
study [Ref 2]. This study concluded that as the liquid reorients and 
accelerates, the free surface begins to assume a curved interface. 
Analytically breaking the liquid into finite sections, (Fig. 3-1) 
reveals that each segment exerts a centriptal acceleration force on the 
tank perpendicular to its velocity vector. Figure 3-1 shows that the 
net centriptal acceleration force exerted on the tank is less than that 
obtained by lumping the total liquid mass at its center of mass. The 
graph in Figure 3-1 presents the mass factor versus percent fill volume 
relationship derived in the referenced study. 


The LAMPS model was used to simulate the Phase I drop tower test data. 
Figure 3-2 shows a comparison of the this model and drop tower test 
data for the z-axis force. This comparison is for drop tower test 22, 
a 15% fill volume case. The model produces a large force spike as the 
center of mass moves through the large change in curvature at the top 
of the ET LOX tank. The integrated effect or impulse applied to the 
tank is shown in Figure 3-3. Again, this result is conservative with 
respect to the test data. , , 
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Liquid Effective Mass Justification 




Liquid Effective Mass Comparison 
Fi.jio'e 3-1. Liquid Effeatiue Maas Development 


Application of a reduced mass factor produces a reduction In the levari 
of forces but re.iults in a nonconservative simulation of the impulse 
applied to the tank. These results are shown in Figures 3-4 and 3-^. 
From these observations, it was concluded that a distributed or multi- 
ple mass representation of the liquid motion was required to more 
closely simulate the propellant dynamics. 
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Figure 6-2. 

LAMPS Model Comparison to Divp Touei^ Test 32 Data^ 
Tank Z Axis Force 


3.2 lJiMPS2 MODEL CONCEPTUALIZATION 


The drop Cower and KC-135 sts were performed with biaxial accelera- 
tions applied to the tank, which resulted in planar reorientation. 
Evaluation of a multiple mass approach could therefore be performed 
using a two-dimens iona representation of Che liquid. This representa- 
tion was obtained by passing a plane (Y-Z) through the axis of symmetry 
of Che «.ank. This produced a two-dimensional slice of liquid that was 
divided into triangular regions (finite elements). The thickness, nor- 
mal to Che paper, varies from element Co element but was assumed con- 
stant within the boundaries of a given element. The total liquid vol- 
ume was Che sum of Che volumes of Che triangular regions. The mass 
within each element was collocated to Che vertices or node points of 
each triangle. This resulted in a distribution of the liquid mass 
throughout the liquid volume (Fig. 3-6). 

Dividing the liquid in this fashion results in two generalized coordin- 
ates for each liquid node: Che y and z positio coordinates of each 
node point. The element •• ■le points were caCego. ized according to lo- 
cation in the liqirid vclu.. Type I node points occur on the free sur- 
face, type IT node points lie on the tank boundary, and type III points 
are chose points interior to the liquid. Under normal circumstances, 
the liquid cannot penetrate the tank wall boundary. At the tank 
boundary, the velocity of the liquid normal to the wall must be zero. 
Therefore, nude points in the vicinity of the wall may only have velo- 
cities tangent to the wall. This constrains the type IT node points to 
follow the tank wall. 
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Figure 3-3. 

LAMPS I4cdel Ccmpariaon to Drop Tower Test 22 Data, Impulse 
of Tank 1 Axis Force 


This constraint is similar to that employed in the LAMPS formulation 
[Ref l]. The difference being that the actual tank boundary is used as 
the constraint surface, whereas the LAMPS constraint surface represents 
the locus of center of mass points, assuming the liquid has a flat 
surface. Test observations have shown that, in fact, the surface 
becomes curved as the liquid circulates the tank. Therefore, the 
true center of mass locations would be different from those assumed in 
the LAMPS model. The distributed mass approach circumvents this prob- 
lem by using the tank boundary as the constraint surface. 

The forces interconnecting the nod® points result from the liquid being 
incompressible and capable of developing internal viscous shear 
forces. The liquid incompressibility manifests itself as a set of non- 
holonoroic onstraint conditions applied to the generalized coordin- 
ates. The incompressibility forces, along with the boundary constraint 
forces, are included in the equations of motion by using the Lagrange 
undetermined multiplier method. The viscous shear forces appear as a 
set of dissipative generalized forces applied to the right side of the 
equations of motion. 
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FicuPB 3”4t HI 1. 

LAMPS Model Comparison to Drop fOMtr Test 22 Data, Tank 
Z Axis Faroe, Mass Factor = ^ lb 
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Figure 3-5. 

LAMPS Model Comparison to Drop Tower Test 22 Data, 
Impulse of Tank Z Axis Faroe, Mass Factor » .15 
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The equacions of motion are numerically integrated with respect to time 
to obtain the new position coordinates of the liquid. These new coor- 
dinates are used to update the finite element geometry. The new geome- 
try is employed to update the liquid incompressibility and viscous 
forces. The equations of motion are again evaluated and the procedure 
repeated. This approach produces a time history of the liquid motion 
which may be directly compared to the photographs of the test data. 
Furthermore, the forces the liquid exerts on the tank can be correlated 
with those measured in the drop tower and KC-135 tests. 
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DEVELOPMENT OP THE EQUATIONS OF MOTION 


The primary objective of developing the tvo -dimensional model was to 
demonstrate the applicability of the finite element formulation in 
solving the liquid motion problems. This model was designed to be cor- 
related with the drop totrer and KC-135 test data. In both these data 
the liquid swtion was planar and occurred in a nonrotating frame of re- 
ference. With these considerations in mind the equations of motion for 
the liquid were developed from Lagrange's equations. These equations 
consider generalised forces and employ nonholonomic constraints. This 
classic form can be expressed as. 


dt 




Qu + b 


kJl 


( 3 -?) 


with the auxiliary equations. 



a 



0 


(3-3) 


where 

q ■ kth generalized coordinate (k»l, 2, ...n coordinates); 

K 

L “ Lagranglan (kinetic - potential energy) ; 
b. ■ coefficients of the Lagrange multipliers; 

■ Lagrange multipliers (i “1, 2, ...a constraints); 

“ the generalized force for the kth coordinate. 
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The Lagrangian (L) for the liquid was assumed to consist only of kinet- 
ic energy terms. The velocity, and hence the kinetic energy for the 
ktb liquid mass point, can be written with the use of Figure 3-6. 


+ Rj (nonrotating frame) 


(3-4) 

(3-5) 


where 


q = velocity of the liquid relative to the tank frame; 

K 


= velocity of the tank frame; 
= kth mass point. 


The positions of the liquid relative to the tank frame were selected as 
the generalized coordinates (all test liquid motion was observed rela- 
tive to the moving tank frame). Therefore, evaluating the partial de- 
rivatives of Lagrange's Equation, 



(3-6) 


Hence, Lagrange's equations for the liquid can be written as, 


» .. (?-7) 

q, = -M, + Q, + b, „ X, 

k K T k kJl i 

This expression, along \ ith equation (3-3) can be recast into a conven- 
ient matrix format, 

[M] {q} . -[M] {R^} + {Q} + [b]^ It] {X} (3-8) 

(nxn) (nxJ) (nxn) (nxl) (nxl) (nxm) (mxm) (mxl) 

[e] [b] {q} = (0) (3-9) 

where 

[e] = diagonal matrix of I's and O's; 

= 1, turns i th constraint on; 


^ii ~ turns i th constraint off. 
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Combining the [e] matrix and the [b] matrix, a short hand notation 
can be employed, 

[b^] = [e] [b] (3-10) 

Before proceeding vith the detailed development of the generalized 
forces and constraint conditions, the techniques used to solve this set 
of equations will be discussed. 

Multiplying equation (3-8) by lM]“l and employing the short hand no- 
tation , 

{q} =- - {R^} + (b^]^ U> + [M]"^ {0} (3_11) 

Differentiating (3-9), 

[bj {q} + [bj |q| = {0} (3-12) 


Substituting equation (3-11) into equation (3-12), {X} can be solved for, 

(3-13) 


{A} 


[bj [M]'^ [bj^ 


[b J {r„> - [b 1 (ql - [b J {Q} 


The values for |x| can he substituted back into equation (3-11) to ob- 
tain the solution for the generalized accelerations. 

3.3.1 Conservation of Momentum 


The boundary nodes are permitted to leave the constraint surface, trav- 
el interior to the tank, and recontact the constraint surface at anoth- 
er location. Likewise, free surface points may contact the tank con- 
straint surface. During these impacts, the momentum of the liquid mass 
normal to the constraint surface is imparted to the tank. Therefore, 
by conservation of momentum (assuming a coefficient of restitution 
equal to zero: the liquid does not rebound), the impulse-momentum 

equations in terms of the constraint conditions can he derived. 




{f} = (Ml 




(3-14) 


3-9 



where 
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f = the impulsive force; 
q = velocities prior to contact; 

= velocities after contact. 

Hence, 

[M] {q}^ = [M] {q}_ + {f} ( 3 .I 5 ) 

{q}_j_ = {q}_ + [M] ^ [b^]^ {f} (3-16) 

Recalling Chat, after contact, Che liquid velocity normal Co Che wall 
is yero. 


1 


[bj {q}_^ = 0 (3-17) 

Substituting equation (3-16) into equation (3-17) and solving for the 
impulsive force, 


{f} = 


[bj [M]"^ [bj^ 


-1 


[bj {q}. 


(3-18) 


Substituting equation (3-18) back into equation (3-16) the final ex- 
pression for the velocities after contact is, 


(3-19) 


{q}+ = {q}_ - [M]'^ [bj^ 


[bj [M]'^ 


-1 


[bj {q}. 


3<3.2 The Definition of the Conditions of Constraint 


Two constraint conditions are imposed upon the equations of motion for 
Che liquid. The liquid, which is assumed to be incompressible, re- 
quires the element volumes to remain constant. This is a time invari- 
ant constraint condition. The number of incompressible constraints 
equals Che total m;- '>er of liquid elements. The second constraint con- 
dition requires the liquid to remain within the boundary of the tank 
(i.e., the tank wall is impermeable). This is a time varying con- 
straint condition, in chat, node points may leave and recontact the 
tank wall. The number of potential impermeable constraints equals Che 
sum of boundary nodes (type II) and free surface nodes (type I). These 
conditions are expressed in Che elements of Che [b] matrix. Although 
Che number of actual constraint conditions imposed may vary with time, 

3-10 


k 


i 
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the size of the [b1 matrix remains constant due to the definition of 
the [e] matrix (equation (3-10)]. This approach greatly simplifiea the 
progranming of the constraint equations. 


Impermeability Constraints - Those node points that lie on the tank 
wall must satisfy the scalar equation, 


n . 
J 



(V20) 


where 

n. = the outward normal vector to the constraint surface at the jth node 
point location; 

= velocity vector of the jth boundary node point. 

The normal vector is defined as the gradient of the constraint surface, 

1 ^, evaluated at the jth boundary node point location. 




z . 
J 



z . 
J 


where V is the differential operator. 


(3-21 


The constraint surface selected to represent the tank wall is a dis- 
torted ellipsoid. Referring to the tank sketch (Fig. 3-7), the con- 
straint surface can be defined as, 

, / \2 O 

ifi = (y/bo)"^ + (z/(ao - aja)] -1 =» 0 

where 

““ ■ “fwd 

■ (“aft - “fwd)/(%ft “f»d) ■ 
bQ = semiminor axis of ellipse (tank radius). 


Using this definition of the constraint surface, the partial deriva- 
tives can be evaluated to obtain an analytical expression for the im- 
permeability constraint of the jth boundary node. 


n . 
J 





z . 

3 Z . 


0 


(3-23) 


3-11 




Figure 3-7. Constraint Surface Definition 
In matrix form. 


bo2 


2ao Zj_ 


(ao - aj Zj) 3 


1 

i’-l 

. 1 

[^4 


^-<o> 


(3-24) 


Equation (3-12) requires 


fKi 3b ^ 3b 3z 
3? 3F 


Evaluating the partial derivatives. 


Hence , 
(3-25) 


3y bo^ 

3b 2^0 + ^ao ai z 

» — ^ - .... - ■ I ■ J 

3 z |ao - ai Zjp 


(3-26a) 


(3-26b) 


and writing in matrix form for the j th boundary point, 



- 

/ 2 \ 

” 

fh! = 

l~\ y i 1 

2ao + 4ao aj z^ 

1. 

IMj 

\bo2j 1 1 

(ao - a, Zj)- ^ 

rj 


(3-27) 


3-12 


i 




r 

CF k\ 




: / 


Constraint Surface Definition for the Slosh Baffles Rpgion - The ET 
tank contains slosh baffles that perturb the liquid as it circulates 
the tank. The concept for incorporating the baffles involves superpo- 
sitioning the distorted ellipse with a (1 - cos) function to represent 
the baffles region. The surface equation development begins with the 
definition of the baffles function. Figure 3-8 shows a schematic of 
the tank baffles and the definition of the (1 - cos) function, rj . 

where 

n = (1-cos -—) h/2 

c = z - + L/2 

h = Baffle Height 

(or Effective Height) 



(3-2d) 

(3-29) 


N = Number of Baffler 




Figure 3-8. Constraint Surface Definition in Baffles Region 


Combining equation (3-22) and baffle function in the region of the baf- 
fles yields the following surface equation. 



(j) = 4,^ + 


(3-30) 


3-13 
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V4> for the baffles region can be written as, 



Defining, 



(3-31) 


(3-32) 


and performing the Indicated partial differentiations. 


9i|) 

e 

3y 


2)^ 

bo 




e 


3z 


3a 

3z 


(3-33a) 

(3-33b) 



0 (4i is a function of z only) 

b 


(3-33c) 




2_ 3n 
bg 3z 


2 3n 

bg 3 


(3-33d) 


where 


, -2ag z. 

3a ^ ^ 3 

3z 


(ag - ai Zj j 3 


and . Ih - /2 tt 


3z 


sin 


In matrix form. 


; 

_ 3a n 

r 3a -Js 

2 3n' 

^ 2 

3n 1, 1 ' 


1^1 

3z bg 

[sz 

bg 3zJ 

bo 

9r“‘‘j 

[M 


(3-34) 


= [b]j {0} 


In a fashion similar to the development of equation (3-27), the time 
derivative of the [b] matrix for the baffles region can be written as, 


3-14 



i 



where 


S^j7.^;vy 


“2^0 ap - 2a 1 z^' 



Propellant /Tank Interaction Forces - The constraint force, X , is 
defined positiv*' (compression) V »n acting inward toward the liquid and 
negative (tension) wuen acting outward. Therefore, when the constraint 
force is acting to prevent wall penetration, X is positive. The units 
of [b] in equation (3-9) are (1/cm), and the units of [n] in equation 
(3-12) are (l/cm-s). Hence, the units of the impermeability lambda's 
in equation (3-13) are. 


ri_ 1_ [-4-] = SbLh = dyne-cm 

Lem gr cmj L s2 J S‘ 


(3-36) 


The units should be dynes. The inconsistency occurs because the defi- 
nition of the impermeability constraint, equation (3-20), used tie nor- 
mal vector, n j , instead of the unit normal, n/lnl. Therefore, the 
impermeability lambdas must be multiplied by the magnitude of the nor- 
mal to obtain the proper units. These forces are resolv<^-i into the 
tank coordinate system by using components of the normal vector. 
Therefore, the total liquid forces acting on the tank are. 


F 

y 


j 





F = 
z 




M = 

X 


j 


2 ( 




(3-37) 


(3-38) 


(3-39) 


Boundary C'^nstraint Relaxation and Enforcement - Because the impermea- 
bility constraint varies in time, criteria for switching the con- 
straint are required. Criteria for enforcement and relaxation of the 
constraint must satisfy the physics of the problem and account for the 
limitations of the numerical integration method. 


3-15 
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The obvious physical condition for enforcing the constraint is that the 
node point reaches the boundary. Substituting the node points position 
coordinates into the surface equation, 

‘‘‘(it- di z) 


(y/bo) 


will yield a value for ^ identically equal to one when the node is on 
the boundary. Coordinates of a point insiJe the boundary will have a 
value less than one and points outside the boundary will have a value 
greater than one. Therefore, by evaluating a criteria for 
enfor-:"‘nent can be established. 


Nunerical integration is not a continuous solution but rather a solu- 
tion at sotae discrete interval of time. T.<erefore, node points ap- 
proaching the boundary (possible enforcement) wil* j'nerally not fall 
exactly on the boundary. This problem can be solved by defining two 
parallel surfaces interior to the actual boundary. Referring to Figure 
3-9, define a* roost interior surface, and ^2 the surface 
closest to the actual boundary. These surfaces are located epsilon 
distances from the tank boundary and £ 2 ) such that, ■ 0.98 
and ^2 ■ 0.99999. A node point (point A, F’^ure 3~9) tnt'jring the 
region between (fij and the cank boundary shall be considered a candi- 
date for constraint enforcement. Penetration of the boundary by 
P'.int A is prevented during the next integration time interval, as fol- 
lows. If the velocity vector of point A, (V^), is directed towards 
the boundary, penetration is likely. Therefore, If tlie angle, 0 , be- 
tween the velocity vef"' 'v and the tank surface normal is less then 
91°, the constraint iw enforced. A node point (point B) entering the 
region between (jjj and the boundary is considered to be in fhe proxi- 
mity of the boundary resulting in constraint enforcement. 

Test observations have shown that the liquid will leave the tank wall 
when encountering the baffles. To account for this phenomenon, the 
constraint forces are relaxed (removed) if three criteria are met. Re- 
ferring to Figure 3-10, these criteria car be defined as, 

1) A<.0, tensile force on liquid; 

2) FT > GF * Fa, the tank/liquid interaction force component in the 
direction of the D'Alembert force. Fa, is grea'er than a factor, 

GF, times the D'Alembert force. Note, the force on the tank la 
opposite the liquid constraint force ; 

3) The node point to be relaxed is adjacent to a free surface node. 

The first criterion indicates the liquid is "pulling" on the tank wall 
and is no longer trying to penetrate. This force should equ«:; or ex- 
ceed the projection of the D'Alembert forces onto the surfacvr r.^nnal. 
The coefficient, GF, provides a mechanism to account for aJhtfive or 
other empirical forces th?. liquid might exert on the tank lgur.,.>rally, 

GF was equal to 1.0 in the correlation analysis). The third criterirn 
prevents a bubble from forming between liquid and the wall. 
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Figure 3-9. Constraint E ■'foraement Criteria 

Incompressibility Constraints - The incompressibility constraint is 
maintained at all times, therefore it is an invariant constraint condi- 
tion. This constraint requires the time rate of change of each finitf» 
elv„ment volume to be zero. The finite element employed in this modnl 
is a constant thickness triangular element that has one node point at 
.ach vertex. Figure 3-11 depicts a typical element with a vertex niim- 
oering of i, j, k.Tf two vectors, a and b, are defined where a is a 
vector from the vertex i to k and b is a vector from vertex i to j , tim 
element volume can be written, 


3-17 


ORIGINAL PACE ?C 
OF POOR QUALII V 



Figure 3-10. Constraint Relaxation Criteria 
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Figure 3-11. Triangular Element Definition 
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■¥*A t"v(axb) = constant 
e e 2 

where 

¥ > element volume (constant) ; 
e 

. = element arv:a; 

e 

t = element thickness 

• ' - i-iM' ^ 

b - (yj - yj|j + (bj - b.|k 

Expanding the vector cross product and dividing through by 1/2 t, 
equation (3-40) becomes, 

(\ - >'j) rj - 'i) - (\ - "i) (>'j - ^'i) ■ 


The time derivative of this expression must equal zero. Therefore 
differentiating with respect to time, 


(‘k - ‘j) * ('i - + (^j - + (yj - \)*i 

^ * (''i ■ 'j) ^ ■ “ 

Casting into the matrix form: [b] |q} , 


ryi>i 








z . 
J 




(3-40) 


(3-41) 


(3-42) 


(3-43) 
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The cime derivative of [b] for the incompressibility constraint is. 


[bl 


[‘k ■ T j - 




I 

- y ' 


I 

- i I 


■*.] 


(3-44) 


These constraint conditions (number of elements in size) were combined 
with the impermeability constraints (number of exterior node points in 
size) to complete the constraint definition. 


3.3.3 Viscous Dissipative Forces 


The generalized forces, Q|(, applied to the liquid are composed of 
viscous dissipative forces. To arrive at an expression for these gen- 
eralized forces, the viscous shear forces within an element were equat- 
ed to the forces applied at the elesient's node points. Consider a dif- 
ferential volume of liquid within the triangular finite element (Pig. 
3 - 12 ). 



where 

y » Shear Strain Rate 
T =■ Shear Stress* 
dx t, the Element Thickness 
= Shear Force 

*In a fluid, devlatorlc shear stress cannot occur without motion; 
shear Is a function of Velocity Instead of displacement. 

Figure 3-12. Liquid Shear Stress Definition 


Tlie virtual velocities of the liquid can be^ defined analogous to vir- 
tual displacements. Therefore, expressing Ydz as the shearing velo- 
city, the virtual velocities of the liquid in shear are, dz. When 
this was combined with the shear force expression shown in Figure 3-12, 
the virtual power (time derivative of virtual work, 6W) of •’’■e dissipa- 
tive forces was obtained. Thus, 


iRefer to Appendix A tor a detailed discussion of the virtual pov^r 
concept. 
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6W ■ -t6y dxdydt 

The virtual power of the nodal point forces must equal the virtual 
power of the internal shear forces. Defining 6 q^as the virtual vel- 
ocities of the generalized coordinates, the virtual power of the nodal 
point forces can be written as. 


6W - F • 6q ■ - t6y dxdydz 
V e 


(3-46) 


where 

F = nodal point forces due to viscous shearing. 

V 

Solving for the nodal point forces, 

6 Y 

F ■ - T TT dxdydz 

V 6q^ 


(3-47) 


The total nodal point forces for an element are obtained by integrating 
over the element volume. 




dxdydz 



(3-48) 


Now, the shear strain rate can be expressed in terms of the liquid vel 
ocities. Thus, for two-dimensional flow. 


(3-49) 

. 3w . 3v 

' 3y 3z 

where 

V = velocity in y direction; 
w = velocity in z direction. 
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PurChennore, from fluid mechanics, the shear stress can be expressed as, 
T * UY, for nonparallel laminar flow (3-50a) 

or 

T = (u + pe)Y, for turbulent flow (3-50b) 

where 

M = liquid dynamic viscosity; 
p = liquid density; 

e = liquid eddy viscosity (not a fluid property) . 

By employing the laminar flow definition of the shear stress and using 
equation (3-49) the nodal point force definition becomes, 

(3-51) 



The expression for the nodal point forces requires definition of the 
velocity field within the finite element. Displacements and velocities 
of the element can be approximated by linear shape functions. These 
shape functions relate displacement and velocities of any point within 
the element's boundaries to the displacements and velocities of the 
element's node points. The area coordinates of the triangular element- 
al region were selected as the approximating shape functions (these are 
linear shape functions having constant derivatives).* Thus, the dis- 
placement field within the element can be written as, 


y =■ N^y^ + N^yj + (3-52a) 

z = N^z^ + NjZj + (3-52b) 

with the condrtion, 

Nf + Nj + \ - 1 (3-52c) 

The shape functions for each node point are. 


*An excellent discussion of the use of area coordinates in finite ele- 
ments is presented in Reference (131. 
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••i-k 

''j'k ■ Vi * [‘i • 'k] 

y + 

1 

[’'k - ^i] 

1 

"3 -k 

[ Vi - 'k''! ^ fk - *l] 

y + 

[>'i - '’k 

11 

\ “ 2A 

[’'I'i - >'i^ ^ 

jy + 

["j ' 

1-1 


where 

A * area of the element; 

V .z = the coordinates of the element's nth vertex, 
•'n* n 


(3-53a) 

(3-53b) 

(3-53c) 


The area coordinates are functions of geometry only, 
velocity field can be written as. 


V . N^y^ + Njyj + Nj^yj^ 


w - + NjZj + 


Therefore, the 


(3-54a) 

(3-54b) 


The strain rates can now be determined, 


3w 

3y 



^ z + — 

3y i 



z + ^ 
i ^ 3N, 



(3-55a) 


3v 

3z 


or 











(3-55b) 


ly ■ 1/“ [ijPj - - ^i) + - 2j)] (3-S6a) 


For convenience, define a local element coordinate system, that origin- 
ates at node i. Therefore, the j and k nodal coordinates relative to 
node i can be expressed. 
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y = y - (3-57a) 

i = 2 - (3-57b) 

By Using this definition of the coordinates, equations (3-56a) and 
(3-56b) become. 


|2 - 1/2A [ 

- \] * 'i(^) - %(^)] 

(3-58a) 

ll ■ 


(3-58b) 


In addition to the strain rates we require 



1 9w ^ 
\3y 3z 


When performing the indicated partial differentiation the rate of 
change of the shear strain rate is determined with respect to the nodal 
point generalized velocities. 




(3-59a) 

3 ^ 

-y,.^A 

(3-59b) 

dy. «v 

3 

y./2A 

(3-59c) 



3 

/- - \ / 

(3-59d) 


b - 

_3_ _ 


(3-59e) 

V <« « CW 

3 ^ 

3 

1 

N 1 
> 

(3-59f) 


By substituting equations (3-58) and (3-59) to equation (3-51) and in- 
tegrating over the element volume, the nodal point dissipative forces 
are obtained. 
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1 
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> 


Fvy^ - 

ty : 

\ / 
'k y** 


Fykj - 

-ty 

\ y/aa 





Fvy, 


FvZj^ = ty z 


tu -y^AA 

j ''h 


where 

Y 


^ -\)* + *'i(^ - >^j) - Vj] 


(3-60a) 

(3-60b) 

(3-60c) 

(3-60d) 

(3-60e) 

(3-60f) 

(3-61) 


The nodal point forces (Eq. 3-60) are the components of the generalized 
force vector, Q , in equation (3-7). 

(V 
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3.3.4 Blement Mass Matrices 


Two techniques were employed to distribute the liquid mass to the ele- 
ment node points. The first technique was a lumped mass approach. The 
mass of each element was calculated (based on the element volume) and 
divided equally among the three node points. Blement volume was deter- 
mined from the element area and thickness normal to the Y-Z plane. The 
thickness was permitted to vary from element to element dependent upon 
the tank width (X direction) at a given element's centroid. This re- 
sulted in elements at the tank boundary having less nuiss than those at 
the center of the tank. Ihis method of mass distribution for an ele- 
ment produces a (6x6) diagonal mass matrix. The element mass matrix 
[equation (3-62)1, has 1/3 of the element's mass lumped in the Y and Z 
direction at each element node point. 


M = fitA 
E 


1/3 


1/3 


1/3 


1/3 


1/3 


1/3 


(3-62) 


The second approach developed a consistent mass matrix for each ele- 
ment. This form of mass matrix was derived from the kinetic energy ex- 
pression for the liquid. Kinetic energy can be expressed by using gen- 
eralized coordinates, q, and the liquid density, P. 


K.E. 



(3-63) 


It is assumed that the displacements of these generalized coordinates 
can he approximated over a finite element region of the volume by a set 
of linear approximating functions. Thus, for one element. 


|q[= [N] {q^}: (3-64) 

where 


(N) - The shape functions over an elemental region [see 

equation (3-53)]; 

|qe{= The nodal point displacements for the element. 

Tlie kinetic energy expression requires the time derivative of the gen- 
eralized coordinates. Furthermore, this representation of the veloci- 
ties should be consistent with the representation of the displacements; 
therefore, recalling equation (3-54) the generalized velocities can he 
written as. 
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Substituting the expression for the generslized velocities into the 
kinetic energy expression and summing over the number of elements used 
to represent the liquid volume results in, 
l^numel 

i >P I L 

dV. 


K.E. - 1/2 


E f ^ 

1-1 




(3-66) 


The element mass matrix is defined from this expression as, 


"• 7 . 


[N]" p[N]dV^ 


(3-67) 


If the expressions for the shape functions [equations (3-53)] are 
substituted Into the Integral, assuming constant density and con- 
stant thickness over the area of the element, the results are: 


Mg * PtJ A 


c/a 




N,N, 
J i 


NjNi 




NfNj 


0 






N^N, 0 
1 k 


0 


Nl\ 




N.N^ 
J j 




0 N, N, 0 
k k 


N, 0 N, 0 N, N, 

k 1 k j k k 


dA 


Now it can be shown that, 


(3-68) 


ff ■ { 


1/12 area, r # i 
1/6 area, r - s 


(3-69) 


Tlierefore, the consistent mass matrix for an element i , 
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1/6 


0 1/12 0 1/12 0 


M - pt A 
e 


0 1/6 0 

1/12 0 1/6 

0 1/12 0 

1/12 0 1/12 
0 1/12 0 


1/12 0 1/12 

0 1/12 0 

1/6 0 1/12 

0 1/6 0 

1/12 0 1/6 


(3-70) 


Note that the consistent element mass matrix produces inertial coupling 
between the element node points. For example, the first row of equa- 
tion (3-70) couples the y coordinat s for the element's three node 
points. Thus, the inertial forces in the y direction at node one re- 
flect the y acceleration at nodes two and three. The total mass at 
each node is decermined by summing the consistent mass matrices over 
all the elements. 


This discussion concludes the development of the two-dimensional equa- 
tions of motion for the liquid. These equations were incorporated into 
a FORTRAN IV program called LAMPS2. A user's guide for the program is 
delineated in Appendix B. The next section of this chapter discusses 
the implementation of the computer prop am and some special numerical 
techniques used in solving the equations of motion. 


).4 COMPUTER PROGRAM IMPLEMENTATION AND NUMERICAL TECHT .QUES 


Several numerical techniques were required to implement the equations 
of motion into a workable computer program. These techniques included 
methods to generate the finite element mesh, initialize the position of 
the liquid in the tank and integrate the equation of motion with a var- 
iable time step. The following paragraphs delineate these techniques. 

3.4.1 Finite Element Mash Generation (MESH2 Program) 


A technique to divide the liauid volume into a set of t> iangular finite 
elements was required. A two dimensional mesh generation program, 
MESH2, was developed to meet this requirement. This mosh program ge- 
nerates a triangular finite element mesh for the axisymmetric orienta- 
tion of the liquid. Asymmetric mesh arrangements were obtained by the 
technique discussed in Section 3.4.2. 
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The axisymmetric orientation of the liquid (Fig. 3-6) is bounded bv the 
tank wall and liquid free surface. The tank wall is defined by the 
constraint surface equation. The free surface location of the liquid 
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can be determined from the tank geometry, fill volume and the accelera- 
tion vector direction. An iterative method to calculate the free a r- 
face location had previously been developed [Ref 1], and was available 
In a subroutine called FLUDC6. The MESH2 program used this subroutine 
to calculate the free surface location of the liquid In the tank. 

Having analytically defined all the boundaries of the liquid, the prob- 
lem is to divide the liquid in a uniform fashion. Due to the curved 
boundaries of the tank wall, this is difficult to accomplish in the 
tank rectangular coordinate system. The problem can be solved by using 
curvilinear coordinates to map square meshes in the plane into 
quadralaterals with parabolic sides in the Y-Z plane. An example of 
this is shown in Figure 3-13. 



5-n Plane Y-Z Plane 


Figure S-IS. Curvilinear Mapping 


This mapping indicates that for every point in the C-n plane, there is 
a corresponding point in the Y-Z plane. Therefore, the plane can 
be divided uniformly to produce uniform divisions in the Y-Z plane. 

The resultant mesh in the Y-Z plane reflects the curved boundaries. 


This method of mesh generation requires the definitions of the mapping 
relationships from the plane to the Y-Z plane. 


y ■ y(5tn) 


(3-71a) 


z “ z(tfn) 


(3-71b) 


The MESH2 program employs a set of "serendipity" shape functions 
fRot 13l derivcii from an eight-nodec* rectangular elemei.t in 
the 5-q plane. In terms of shape functions, the mapping rt let ioosh i ps 
become , 
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*(S.n) “ 2 zcj, 

1 

(3-72a) (3-72b) 


where 


n 





Nl 

- +1.0 

(1 - 0(1 - n)(l - 25 - 






N2 

- +4.0 

(5)(1 - 5)(1 - n) 

7 

6 

s' 



N3 

- -1.0 

(5)(l - n)(l - 25 + 2n) 

( 1 A 


A 1 



N4 

- +4.0 

(5)(n)(l - n) 

’ ' o 


H ' 



Ns 

- -1.0 

(5)(n)(3 - 25 - 2n) 

,JL- 

2 

3 



Ne 

- +4.0 

(5)(1 - n)(n) 


-1.0 (1 - 0(n)(l + 2C - 2n/ 
+A.0 (1 - 0(n)(l - T„ 


Thp (yc^,zc^) coordinates arc eight points .hpt define the boundary 
of the region toat is being napped by tne {,-n plane. The three points 
along any side are used to define the parabolic curve in the Y-7. plane 
used for that side. The mesh definition is obtained by defining the 
number of divisions along aides 1 and Ttie plane in then divided 
into the specific number of divisions and the £-ti coordinates are 
obtained. The ?-n coordinates are then transformed, by using the 
ahapo functions, into t*’e Y-Z plane [equation (3-72)1 

Figure 3-14 shews two arran 'cements f the boundary nodes used to define 
two raesh orientations. In either case, the node point-side definition 
is, 


Side I - Nodes 1.2,3 

Side 2 - Nodes 3,4,5 

Side 3 - Nodes 5,6,7 

Side 4 - Nodes 7,8,1 

After mapping, if the mesh boundary points do not match the boumla.'y 
location as defined by the distorted ell’ ’oidal constraint surface 
(due to parabolic mapping), the boundary points are moved to mate the 
constraint surface defi.iition. 

The triangular elements are formed by selecting the shortest d'agc.nal 
in each square after mapping to the Y-2 pi ?n . If a saries of rectan- 
gular elements are erected with equal diagonals, the direction of the 
diagonal is alternated from one element to the next. 
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Examples of the type I and type II mesh arrangements are shown for a 
lOX full ’■-ank in Figures 3-15 and 3-16, respectively. Both mesh types 
were evaluated in the correlation analysis. 
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Arrangement for Meah Generation 


3.4.2 Liquid Initial Position 

The mesh generation program considered only the axisymmetric orienta- 
tion of the liquid. Initial positions of the liquid, other than the 
axisymmetric case were required for the correlation analysis. These 
asymmetric orientations were determined by an iterative scheme that 
minimized the potential energy of the liquid. 

Figure 3-17 shows a sketch of the tank tilted with respect to the grav- 
ity vector. He'-ce, the axisymmetric orientation is not the equilibrium 
position. There is a delta potential energy (4PE) between the center 
of mass position of the axisymmetric orientation and the true equilib- 
rium position. Minimizing this APE yields an equilibrium state for the 
liquid in the tilted tank. 

Calculation of the potential energy requires a reference datum. This 
datum will vary with the orientation of the tank with respect to the 
initial acceleration vector. This makes it difficult to calculate the 
potential energy of the liquid. An alternative approach is to calcu- 
late the kinetic energy of the liquid as it moves from the axisymmetric 
orientation to the equilibrium position. The kinetic energy peaks as 
minimum potential energy is achieved. Furthermore, if the liquid ve- 
locity is killed at the point at which the kinetic energy is a maximum, 
the liquid remains at rest, since it is at a minimum potential energy 
state. 
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Figure 3-1?. Aeyrmetrio Liquid Orientation 


Tlie method by which Che kinetic energy is determined relies on numeri- 
cal integration of the equations of motion given the initial accelera- 
tion vector. The direction of the acceleration vector is defined by 
specifying a tilt angle of the tank (Y) with respect to the gravity 
vector. This is shown in Figure 3-17. 


The kinetic energy tisM history produced by this technique is not con- 
tinuous. Therefore, the point at which the slope of Che kinetic energy 
time history is cero (a minimum or suiximum) can not be determined ex- 
actly. The rate of change of KB can be expressed by a difference equa- 
tion based on Che current kinetic energy (KBc) and the kinetic energy 
at Che previous integration step (KE|^)> 


^ (3-73) 

dt - At 


where: At is Che integration tisw step. 

The LAHP82 program requires that this approximation of the slope of the 
kinetic energy timt history be less Chan a small nua^er, cke, for 
convergence to an equilibrium state* The program facilitates conver- 
gence by killing the liquid velocity as the kinetic energy passes 
through a maximum. Figure 3-18 shows an example of this stratagem. 
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Figure 3-18. 

A Stratagem for Liquid Position Initialisation Based on 
Kinetic Energy 
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The cop graph ahowa an energy dlagraai of the liquid aa it approachea 
the equilibriua poaicion. The lotier graph depicta the kinetic energy 
tiiae history. 

An additional capability included in the convergence procedure is the 
ability Co aove the liquid to Che deaired tank tilt, in increiaenta of 
the tilt angle. The angle of tilt ia divided into NGAHMA aegaenta and 
the liquid equilibriua poaitiona are deCerained for each aucceaaive 
aegaent. Tbia capability ia uaeful for large anglea of tank tilt 
(Y > 20®). 

3.4.3 Variable Integration Step Algoritha 

The purpoae of Che variable step integration algoritha ia to prevent 
tank wall penetration by a free aurface or boundary node point. This 
penetration can result froa too large an integration interval. This 
penetration can be prevented by reducing the aagnitude of the integra- 
tion step sise. The integration interval required to prevent penetra- 
tion, DTNEH, can be calculated by estimating the new positions of the 
nodes based on their current state definitions. 

Referring to Figure 3-19, the estimated position, Rg, of a point may 
be obtained through Che following kinematic relationship. 
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R^j - Rp ♦ Vp DTNKW ♦ 1/2 Ap (DTNEW)* (5-74) 

Where: 

R * Current Poeiclon of Point pi 
P 

Vp • current Velocity of Point p; 

A “ Current Acceleretlon of Point p. 

P 

The eetimeted position for e given value of At can he euKerituted into 
Che constraint surface equation to evaluate wall penetration, the pro- 
cedure used in the progran starts the estisMte with a mininum value of 
integration step, isinimusi value is obtained by dividing 

the input DTMAX by an integer, N (N is currently set to S in the pro- 
gram' . If the position estimate is inside the boundary, At is incre- 
mented such that, 

At - (At , ) (INCR) ; INCR - I, 2, . . N. (3-7'^) 

nln 

The new value of At is substituted into equation (5-74) and boundary 
penetration is evaluated again. The largest fraction of DTMAX that 
does not cause wall penetration is the mavinmm integration step to be 
used for that node. All free surface and boundary nv>des are evaluated 
in this isanner and the minimum DTNBW of the set is selected as the nest 
integration step site used in the solution of the equa ions of sK>tion. 
This technique has proven verv successful in preventii.g wsll penetra- 
tion. 
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A priaary objaetiv* of this otudy wot to dovolop o toot-worifiod ono- 
lyticol wodol coptblo of tiaulatint IT largo anplituda tloah. Accov 
plithaant of this objactlwa raguirad eorralation of tha Photo I and 
Photo II toot data to tha two-dlaanaional analytical nodal. Thraa eri- 
taria wara aalectad to attatt tha accuracy of tha tinulation: 

1) Correlation of tha liquid notion charactar, with that obtarvad in 
tha high-tpaad noviat (i.a., liquid poaition and fraa turfaca ahapa 
tine hiatory correlation to the teat data); 

2) Force/nonent nagnitude and phaae correlation; 

3) Correlation of the total iaqiulae applied to the tank. 

The firat criterion providea a qualitative eawcaanent of the overall 
liquid notion ainulation. Thia criterion ia avaluatad by co^»aring the 
photographic inagea of the teat data, at apecific tinea, with the con- 
puter-generated graphica of tha ainulation. 

The force/nonent correlation (second criterir i) provides a quantitative 
assessment of analysis accuracy in representing the tank/liquid inter* 
action forces. 

The third criterion provides a measure of model conservatism. The 
total impulse or double integral of the force and mo me nt tine histories 
indicates the degree to which the liquid notion will affect BT rigid 
body dynamics (Ref 4]. A total i^>ulae of the simulated forces greater 
than the teat-derived total impulse indicates the analysis is conserva- 
tive in predicting BT response to liquid notion. The point mass model, 
previously discussed, has been shown to be overly conservative by this 
criterion (Kef 4J. 

Several model parameters were varied in the correlation study. These 
paraiseters include the finite eleamnt geoswtry, analytical representa- 
tion of the baffles region, and the nuanrical integration step sise. 

The finite element geosntry ia defined by the mesh type (I or II), mesh 
sise (number of eleswnta), and the mass matrix formulation (lumped or 
consistent). The analytical representation of the baffles region was 
evaluated by varying baffle height, spacing, and number of baffles rep- 
resented. Analytical convergence to the test results is a function of 
the numerical integration step sise. The solution sensitivity to thia 
parameter was assessed by varying the integration step site for a given 
model geometry. 

The analytical model was correlated with the Phase 1 drop tower test 
data and the Phase II, KC-135 test data. 
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Initial drop towar taat/analytical corralationa indicatad an offaat ba** 
twaan the analytical predictions of the Y-axis force and teat Maaura- 
■ents. Figure 3-20 presents a coe^rison of the test and analytically 
derived forces for test 22. The drop toner test forces reflect the 
acceleration transient frosi the 1-g environaent to the applied test 
accelerations. The test force initially approaches the analytical re- 
sult, then trails off to a value less than the analytical data. 

The force traces in Figure 3-20 consist of longitudinal and lateral 
components of the test forces resolved to the tank Y axis. Investiga- 
tions of these components have indicated that the difference in the 
test and analytical force profile is a result of experimental error in 
the lateral force measurement. Figure 3-21 shows a copy of the strip 
chart measurements for test 22. The lateral force measurement (F^) 
exhibits drift before drop capsule release and erratic behavior for the 
first 0.15 seconds of test time. Study of the longitudinal forces and 
photographed liquid isotion indicates that the suspect lateral force 
behavior is not a function of the propellant dynamics. It is hypothe- 
sized the lateral force experimental errors are a result of the high- 
sensitivity setting of the lateral force transducer and test fixture 
construction. 

Based on the above observations, the drop tower test data were adjusted 
to match the analytical data at a point after the transient event and 
before any significant propellant miction. This adjustsient was accomp- 
lished by the following steps: 

1) The values of analytical and test forces and moments were deter- 
mined at t « 0.15 seconds; 

2) These values were subtracted such that all forces and moments (test 
and analysis) were zero at t > 0.15 s to facilitate correlation of 
liquid motion Induced forces; 

3) All force correlations and Impulse calculations were initiated at 
t - 0.15. 

Analytical correlation with the drop tower data was adversely affected by 
the technique used to determine test accelerations. Drop tower accel- 
erations were calculated from high speed movies as described in Refer- 
ence 3. Consequently, in the drop tower test /analytical correlations, 
which follow, the analytically predicted liquid motion lags the test 
data slightly. However, analytical/test phasing correlation for the 
KC-135 tests was excellent. Accelerations were directly measured during 
KC-135 testing. 

Analytical correlation to the KC-135 aircraft te>.. data was initiated 
at the time the aircraft Z acceleration passed through +1 g [Ref 4]. 

In several Instances, force transducer data channels were saturated for 
a period of time after the +l-g mark. In these Instances, the correla- 
tion analysis was initiated after the data saturation period. No addi- 
tional adjustments were required, since actual accelerations were measured 
during KC-135 testing and the use of strain gage load cells (and larger 
test forces) resulted in more accurate force measurements. 
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Figure 2-20, Teat/Analytical Correlation: Tank Y-Axia 

Foroet Drop Tcnoer Teat 22 


3.5.1 Evaluation of Model Parameters 


Two mesh arrangements were evaluated in the correlation analysis. Both 
lumped and consistent mass matrices were used in the two arrangements. 
Figure 3-22 shows a liquid motion profile comparison for the two mesh 
arrangements using a lumped mass matrix (drop tower test 21 conditions) . 
The liquid motion comparisons for the consistent mass matrix are shown 
in Figure 3-23. 

The lumped mass, mesh type I analysis produced the isost uniform flow of 
the liquid. The second best liquid motion profile was produced by the 
type II, lumped mass formulation. The consistent mass formulation, 
which couples the inertial forces within the liquid, produced the poor- 
est simulation of the liquid motion. Results of these four cases for 
drop tower test 21 conditions are shown in Figure 3-24. Note that good 
liquid motion correlation results in good force correlation. 

It is postulated that the poor simulation quality for the consistent 
mass matrix is a result of the inertial coupling in the equations of 
motion. This inertial coupling appears to require a smaller integra- 
tion step to achieve solution convergence. Since the computational 
costs for the consistent mass formulation exceeded thoae of the luaiped 
mass formulation (for the saaw integration step sise) and correlation 
was not improved, this method was not pursued further. 

A numerical integration convergence atudy was performed using the 
lumped mass formulation. The liquid motion was siimilated using the 
sasie mesh arrangement and sise for various integration step sises 
(DTMAX). Figure 3-25 shows the test/analytical correlations of the 
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tank Z-axia forca and total iapulaa for a KC»13S aircraft taat condl> 
tlon. Thraa integration atap aiaaa vara avaluatad in tha corralation, 

Aa tha intagration atap alta waa raducad, tha aagnituda of tha analytical 
forcaa convargad to tha taat data. Accoapanylng tha aagnltuda reduction 
ia a lengthening of tha forca duration, thua, tha total iaipulaa waa not 
significantly affected. Furtharamra, tha liquid notion profilaa pro- 
ducea for theae intagration atap aiaaa uara vary ainilar. 

% 

Tha amdal traa correlated to ainilar drop tower and RC-135 airera'- taat 
data. Thia correlation produced additional validation of tha taiti 
acaling approach (Chapter 2.0). Tha analytical nodal produced ainilar 
raaulta for tha l/60th- and 1/10-acala taata. Figure 3-26 ahowa forca 
and total inpulae correlation to drop tower taat 5 (lOX fill) and 
KC-13S aircraft teat 1.2.2 (lOZ f'U). In both caaaa, tha nodal corre- 
latea well to the teat data. Liqu.J laotion correlation to tha teat 
photographa ia shown in Figure 3-27 for tha drop tower taat data and 
Figure 3-28 for the KC-13S aircraft teat data. Thia correlation indi- 
catea the model should produce an accurate simulation of full-scale ET 
propellant dynmaics. 

The effect of fill volume on the model was evaluated. The model waa 
used to simulate fill volumes varying from 2Z to ISZ. The saem siae 
mesh and arrangement waa used in all cases. Figure 3-29 shows the 
force and total impulse correlations for SZ and ISZ fill volumes. 

Figure 3-29 presents the correlation for F^, Fy and the moment (N^) 
transferred to the scaled (l/60th scale) location of the external tank's 
CG. This allows for an assessmant of the differences (test/analysis) 
as related to overall external tank dynamics. Reference 4 (chapter IV) 
used the same technique for evaluating slosh simulations for the first 
Shuttle flight (OFT-1). Figure 3-29 shows that the force magnitudes nre 
in general, conservative. However, the Integrated effects (double 
integrals) correlate very well with the test data. 

Figure 3-30 and Figure 3-31 show test /analytical liquid motion comparisons 
for 2Z and 15X fill volumes. The analyses do not correlate cs well at 
the 15Z fill volume. This indicates a finer mesh size should be used 
with higher fill volumes. 

The baffled tank configurations were analysed with the two-dimensional 
model. A multiple baffle approach, which varied analytical baffle 
height, was evaluated. The multiple baffle arrangement consisted of 
eight baffles represented by the superposition of a (1-cos) function 
with the distorted ellipse constraint surface. This approach resulted 
in entrapment of the finite elements in the baffle region and other 
numerical difficulties. Tho result was poor correlation of the model 
with the test data. 

A single-baffle approach, which represented the baffle region as one 
baffle (centtal in region), was evaluated. The analytical height of 
this single baffle equaled the actual baffle height. The liquid motion 
correlation to a baffled KC-135 test is shown in Figure 3-32. The cor- 
responding force and total impulse correlations for the tank Y-axis and 
Z-axis forces are shotm in Figure 3-33. The analytical flow leads the 
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test dst9, but in general correlates veil. This is the preferred 
approach to baffle representation at this time. More work needs to 
be done in this area. 

In conclusion, the ttfo-disMnaional siodel has simulated accurately the 
liquid motion in scale models of the external tank. The degree of 
accuracy can be controlled by the choice of mesh sise and numerical 
Integration Interval. The model produces conservative results that 
are more representative of the liquid motion dynamics than those pro* 
duced by the single*polnt mass analog. 
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Figure 3-23, Mesh Evaluation: Consistent Mass Matrix; Motion 
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Drop Tower Test 21 (5X Fill) 



Drop Tower Test 22 (ISZ Fill) 
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4.0 


mREK“l>TMENSIONAL FINITE ELEMENT MODEI. (IJWPS'^) 


The flnalveis of the .separntion ot the ET from the Orbiter during RTLS 
requires t^w solution of the general three-dimensional equations of u.o- 
r ion for Che propellant. Based on »he correlation results for the two- 
dimensional model, a three-dimensional finite element model was formu- 
lated to analyse the propellant dynamics during RTLS. This model cou- 
ples the ET rigid body dynamics with the propellent dynamics to provide 
a closed-loop analysis of the ET system. This chapter discusses the 
development of the equations of motion for this closed-loop system. 

The ET rigid body dynamics are developed first, fo'.lowed by the equa- 
tions that couple the propellent to Che ET. The coupling equations 
section includes a discussion of the interaction forces between th,» 
liquid and tank. A three-dimensional finite element is derived to mod- 
el the liquid incompressibility and viscous shear forct-S. The numeri- 
cal cechniquea employed Co solve Che equations of motion and the cor- 
relation of the three-dimensional model are included 'i this chapter. 

4.1 EXTERNAL TANK RIGID BODY DYNAMICS 

The rigid body equation of mot’ on fov the ET were developed from 
Lagrange's Equations assuming the Lagrangian only consists of kinetic 
energy terms. These equations included generalized forces and nonhol- 
onomic constraints applied to the ET. 


d__ 

dt 


eu-(gi 


Q. + b 


k« 
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Six generalized coordinates were selected to describe the Ei’ motion. 
These coordinates consist of three displacements of the iIT center of 
mass and three Euler angles, which relate the ET body axes to the iner- 
tial axes. Therefore, the state sector for the ET rigid bodv is, 

' X 


Y 

Z 

ijj 

0 

♦ 


Figure 4-1 shows the ET center of mass relationship to the inertial 
frame. The cenCe ot mass is positioned by the vector Rj. , and the 
orientation of the ET body axes with resoect to the inertial frame is 
defined by the direction cosine matrix [y]- The dirtction cosine ma- 
trix is a function of the Euler angles. Figure 4-2 delineates the 
Euler angle sequence employed in this analysis and the derivation of 
the [y] matrix. 
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Figure 4-1, ET Center of MASS Coordinate Definition 
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Combine 4- 5a, b, and c 



(4-6) 


where 
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ceci^i 

-C(|iSec4»+bii»S4i 

se 

C(JiC0 

-ces\l; 

C(t)S0SijJ+S<J'Ci|' 


S(j)S6Ci[i-l-S4-C<j> 

-S(j)C0 

-S4>S0S4;+Ci})Ci*) 


*r » S * 

Fujure 4-2 Euler Aryjle Sequence 
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Having defined Che aCate vecCor«Che kinetic energy expresaion for Che 
ET can he ohcained. Thia expreaaion cunaiata of the CranalaCional kin- 
etic energy of the center of maaa plus the rotational kinetic energy 
about Che center cf raaaa. 


T - 1/2 {ic Y Z) [Ml < Y >+ 1/2 {u. u) 

X y 


where 

* * • . • 

X, Y, Z * inertial velocities of GT center of mass 

[N] * roasa matrix of KT 



(4-8) 


***x* " bodv axis instantaneous angular velocities 

(nonholonomic velocities) 


(ll ■ ET inertia matrix or 

censor 

inertia tensor is 

defined 

as , 


Ixx 

-Ixy 

-Ixx 

[I] - 

-lyx 

lyy 

-lyx 


-Ixx 

-Ixy 

Ixx 


where 

Ixx - /yP(y^+z^) 
Ixy - oxyd¥ 


mass density 


(4-8a) 


All other terms of this inertia tensor can he chtained hy cyclic permu- 
tation of the indexes x, y, z. 

Substitution of the tirst term of the kinetic energy expression into 
equation (4-1) results in the equation of motion for the translation oL 
the ET center of mass (neglecting constraint forces at this time). 


M 0 

M 

0 mI (4-4) 


The rotational kinetic energy of thr ET most he written in terms of the 
Euler angle gtneriliz d c«>ord inates . ’fhe instantaneous angular veloci- 
ties can he expressed as prelections of the Euler angle rates onto the 
hodv ax*>s. Using Figure 4-2 to ilefine these pro lections. 
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Substitution of equation 4-10 into the second tern of the kinetic ener- 
gy expression yields, 


(4-11) 


- 1/2 U' 0 [ti] (11 lel\ e 



n^e p.irtiul derivatives of this kinetic expression with respect to the 
Euler eagle generalised coodinetes, q j . (j • 4, S, 6), results in, 


. . . 

-- — - 1/2 lii» 0 

■’‘'j 


>'q 


IT 



The inertia tensor is s symmetric matrix, therefore 


^ 

Oq 


li<i e <f) 


.1 



(4-13) 


Define a matrix, [Pj], which contains the partial derivatives of the fS] 
matrix with respect'’ to each Euler angle coordinate. 


[ P . ] ■■ i 1^' 1 1 


'Vi 


, J - 4. 5, 6 


Assembling these matrices into[p], the change in kinetic energy 
with respect to each Euler angle coordinate can be written as. 


(4-14) 


)T, 

4 k 

;-)V 

;)T, 

.Vo‘ 

;)T„ 

d 


(PI m [t<i V' 


i'l 


(4-11) 


4-0 



ORIGINAL PAGu !S 
OF POOR QUAIIT^ 


The partials of the kinetic energy with respect to the general iied 
velocities of the Buler angles arei 



(4-16) 


Stibatitucion of equations (4-15) and (4-16) into equation (4-1) results 
in the homogenous equations of motion for BT rotat:.on, 


(f^r m (Hi 


Irl 
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iKi [II (Hi + inr m (Hi - (PI (11 (Hi 
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T*u> trsns format ions t*»at relate the body axis forces end torques to the 
generalised forces, Q^^, can he obtained through the virtual power 
method. The virtual power of the impressed forces and of the body axis 
velocities is. 



(4-18) 


vdiere 


jll , jv, 

- virtual 
system. 


• virtual 


velocities measured in the body axis 
rotational ratea. 



» fy » ^2 
"'• x » » "^2 
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■ forces measured in the body axis system, 

= torques measured in the body axis system. 


The virtual power in terms of the generalised velocities and forces is , 


5W 


{6X 6Y «Z 66 6i^f<Q3 


(A 19) 


The body axis velocities are related to the generalized velocities 
through the [Y] and the f6] matrices, 


thus , 


{u V w} « {X Y Z> tYl 


(4-20) 


and 


{u) u 

X y z 


(4-21) 


Substitution of equations (4-20) 
suits in, 


6W 


(6X 6Y 6Z 64i 60 6(J} 



and (4-21) 


(cl 



into equation (4-18) re- 


(4-22) 


Comparing equation (4-19) to equation (4-22) the expression for the 
generalized forces is obtained, 
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(4-23) 


These generalized forces were incorporated into the right-hand side of 
equations (4-9) and (4-17) to complete the rigid body equations of mo- 
tion for the BT. 


4.2 PROPELLANT COUPLING TO ET EQUATIONS OF MOTION 


The position of the i th liquid mass point can be specified by using 
the vector relationship shown in Figure 4-3. This states that the ab- 
solute motion of the liquid can be obtained from the sum of the LOX 
tank motion, plus the motion of the liquid relative to the LOX tank. 

The position coordinates of the liquid relative to the tank reference 
axis were selected as the generalized coordinates for the liquid. This 
selection of generalized coordinates results in coupling the ET rigid 
body dynamics to that of Che liquid through the mass matrix of the sys- 
tem. This approach simplifies the expression for the interaction con- 
straint forces between the liquid and tank. Thus, the state vector for 
the liquid and ET is , 

<Sa- V "zl V- ’zn- 


where 

n * the number of node points used to represent the liquid. 
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The Inertia coupling between the liquid and ET rigid body is derived 
from the principal of virtual work. 

/^6Rj(pR^ - f^) d¥ - 0 (4-25) 

where 6R^ •• Virtual displacement of 1th mass point 
R^ * Inertial acceleration of ith mass point 
f^ External forces on ith mass point 
p >■ Mass density 


Solution of equation (4-25) requires 
Figure 4-3, and assuming a rigid ET, 
eratlon of the 1th liquid point are, 

• • • 

R^ - R^ + (I^R^ 4- ii(q^ + q^ 

and ^ ^ • _ ^ 

- R^ + (uXR^ + ujX (uiXR^) + q^^ 

(\) 

where (u)} » < u) >, equation (4-8). 

) 


definition of 6R^ and R^. Using 
the inertial velocity and accel- 

(4-26) 

• • 

+ 2udCq^ + aflCq^ + (4-27) 


The above expression assumed the same basis for all the vectors. In 
general, the tank coordinate system will not be alllgned with the ET 
center of mass triad. Furthermore, the center of mass triad rotates 
with respect to inertial space. The [ y] matrix, equation (4-6), re- 
lates the body axis triad to the inertial triad. The tank triad (T) 
can be related to the body axis triad (C) by. 



where 


(4-28) 


I", J", k" are the tank axes base vectors. 

The R.J, and vectors can be defined in terms of their components. 


and 


Kf - L‘' )' k'j jty I 

■ Li" J" k"J I 5; I 


body axis system 


t k axis system. 


(A-29) 


(4-30) 
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Combining these vectors, the total distance from the ET center of 
mass (empty) to the ith liquid point is. 



h 

-U’ y k'j! 




U' J' k'J j 

h* 

h^ 



r H h' 

^ i 
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i 
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The angular velocity cross products (in terms of the body axis triad) 
can be rewritten in matrix form using this total distance vector. 
Hence, 


* U' j' k’J [r^ 



(4-32) 


where [ r ] . 

s 1 



i 


Furthermore, the centripetal and corlolis accelerations can be written 
in a similar fashion. 


tJC(a)X(R,j, + q^) = L i' J’ k’JCn] 


and 


2a)Xq^ » 


2 Ll’ j’ [^R^] 


where [ fi] 


0 -lu 0) 

z y 

0 -oi 

Z X 

-u (i) 0 

y X 
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(4-34) 


The ET body axis angular velocities, u, and angular accelerations, u, 
can be written as functions of the state variables (ij^, 6, ()>) and 
their time derivatives. Recalling equation (4-10), 


1 

1 • 

1 u 

1 X 

1 

1 U) 

•= [6] 0 

I y 

1 « 

1 ^ 

I 

l z 


(4-35) 
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Differentiating, the angular acceleration of the body triad is- 
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Substituting equations (4-28) thru (4-36) into equations (4-26) and 
(4-27) we obtain the inertial velocities and accelerations of. the 
liquid in terms of the ET and liquid generalized coordinates. 
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The virtual displacements, 6R., are related to the generalized coordi- 
nates through the kinematic relationships of equation (4-37)*. Thus, 
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U ^ Kj 
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^ ' 6^ 


(4-39) 


*See Appendix A, Virtual Power Concept 
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Having obtained axpreaaions for the acceleration and virtual 
displacements of a point within the liquid, the virtual work ex- 
pression can be evaluated to obtain the llquld-ET inertia coupling. 
The virtual work integral (eq. 4-25) Is replaced with a summation 
over n liquid mass points representing the mass of the liquid volume. 
Furthermore, this summation can be separated into acceleration and 
velocity components of the generalised coordinates. Thus, 


Z Lfi5i<Sx6aJ 
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(4-40) 


where liquid mass 

X - X, Y, Z 

a - i|;, e, ^ 


[m^] - mass matrix for the ith liquid mass (3x3) 


The coupled mass matrix is obtained by combining the first summation 
of equation (4-40) with the mass matrices defined in the ET rigid 
body equations (4-9) and (4-17). Thus, 

•• t« •• •• *• 



-n 


a 






[m] “ Sym E [m ]+[m] ' Z [m.][Y]^[r ] [6] 

1-1 ^ ;i-l ^ ® 1 




? C63^[r ]^[.,3[r 1 [B3 

Sym 

Sym 

ll-l “ 1 ‘ " 1 

+ [B3’'[i3[B3 
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I 


4-13 




1 





or In short, 
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The second surjiatlon of equation (4-40) represents the liquid inertia 
forces which are a function of liquid and ET generalised velocities. 
These forces are treated as generalised forces applied to the right 
hand side of the equations of notion and are combined with the ET 
Inertia forces which are a function of the ET generalised velocities 
'4-17). Thus, the coupled equations of motion (neglecting external 
and constraint forces) are. 


P l%] 

k=:% vj (s) <*- 

where, 2 [^R^][£l][^R^][m]j^Hj^ + [^R^,]|^[n]Cr^]^[ 6 ] + [r^]^[ 3 ]j Cn]p 


(4-42) 


LX-J - • r 

2[,R,][!l][<.^3[«]„i„ ie] + [r.] 

- L n n J • 

[Q^] - [xf 

‘ref f £ [t ]■'[■»],[ 2ta][-i>-]<Si> raJCrj [e]{i) + [x ] [e]ill 

CQ„] - r ‘ L ‘ ‘ JJ 

[+ [efCilCe] + [e 3 ’'ti][e] + [P 3 [i 3 ts 3 j{i) 

This completes the development of the inertia coupling between the 
liquid and ET rigid body dynamics. 

4.3 SOLUTION OF THE COUPLED EQUATION OF MOTION 

Before proceeding with the detailed derivations of the constraint and 
liquid shear forces, the general procedure for solving the coupled 
equation will be discussed. This discussion includes the development 
of the equations to conserve the momentum of the coupled system. 

The generalized forces applied to Che liquid and ET, as well as the 
constraint forces, can be added to equation (4-42). 


«HS:"Hx';"sal (-=) \% 

X= XX Xa ) M \ 

„ M M (at Q 

a- av oa V ' a 


1 r T r T " 

!= Qe= *»= 

X ^ 

I Q b 

Ot L . L a J 


(4-43) 
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(4-AA) 


where [Q„-], “ Excemsl end dissipative generalised 

ba EX BQ 4-^ t-hm. m„A 01' 


forces applied to the liquid and ET. 


The first step In the solution process requires that the momentum 
of the system be conserved in accordance with the constraints Im- 
posed upon the system. AnsTtgous to the development presented in 
Chapter 3.0, the Impulse-momentum expression for an Inelastic impact 


can be written as. 


m M 

"sxk- 

1-^1 


CXI 
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TT„ 

az 
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R' 

► - 

SNxiV 

“oH "cio 

m * m 

ll\ 


‘■x 


{f} 

(A-A5) 


where (-) subscript denotes velocities before contact and subscript 
denotes the velocities after Impact. Tl»e Impulsive force at Impect is 
represented by {f}. Furthermore? after the liquid node point has con- 
tacted the wall. 



(A-AO) 


A two step solution process, wWch required inversion of the system 
mass matrix, was employed in Chapter 3.0. The system mass matrix is 
now highly coupled and containa a large number of non-zero terms. 
Therefore, an alternative solution procedure is employed. Equation 
(A-A6) is appended to equation (A-AS) and the impulse forces {f} are 
placed on the left hand side. Thus, 
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The "plus" velocity states are obtained by solving th«je algebraic 
equations with a Gaussian elimination and backward substitution 
procedure. The (e-1) diagonal terms permit decomposition when a 
constraint becomes zero (e«0) . 


Having obtained the generalized velocities of the liquid and ET 
after the liquid Impacts the wall, the velocity dependent forces 
can be evaluated and the coupled equations of motion solved. 
Equations (4-43) and (4-44) are cast into a format similar to 
equation (4-47). 


M„ 


b!l/ 


ZOl 

w 1 




M „ M 
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-Xr-.XX 

X“ 


M „ M 
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b^ 

o“ ax 

oa 

a 1 

b- b 

H '’x 

’’a 

(e-1) \ 
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0 0 0 
0 0 0 
0 0 0 


^5 ‘’x ’’a 


i \ 

X 

m 

a 

(4-48) 


Note, the left hand side coefficient matrix has remained the same, 
therefore, the state solution only requires evaluation of the right 
hand side and a new back substitution. 

The constraint forces between the liquid and tank will be shown in 
section 4.4 to be only functions of the liquid generalized coordinates. 

Therefore, [b^] and [b^] are null matrices. This simplifies the 
expression for the velocity change after impact (equation 4-47) and 
the constraint force expression (equation 4-48). Note, however, the 
Impact of the liquid onto the tank wall still affects the ET velocities 
through the inertial coupling terms. Likewise, the constraint forces, 
A, reflect the ET generalized forces and momentum through the inertial 
coupling. This coupling provides closed-loop feedback between the 
liquid and ET dynamics. 


4.4 THREE-DIMENSIONAL IMPERMEABILITY CONSTRAINTS 

The liquid at the ET wall is constrained to move tangent to the wall. 
The previous section stated that the auxiliary equation, (equation 
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4-39), which defines these constraint conditions, did not directly in- 
volve the ET velocities. This section will develop these constraint 
expressions and define the analytical representation of the constraint 
surface. 


Consider two points separated a distance of A , one point on the tank 
boundary the companion point interior to the boundary. Define the 
boundary point as point A and the interior point as a liquid siass 
point, point B. Referring to Figure 4-4, the absolute positions and 
velocities of both points can be obtained. 


R 

a 


R + 


+ q 

T 


and 


(4-498 ) 


(4-49b) 


R »R +R_ + a)XR„ + q +o)Xq 
a c T T a 


R^ = R^ + R^ + CO X R^ + q^ + a, X q^^ 


(4-50a) 

(4-50b) 
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Figure 4-4, Definition of Three-Dimensional Tank Boundary Constraint 


The velocity of point B, the liquid mass point, relative to wall 
is A or, 


• • • • • __ 

A = R - R = q, + (i! X q, - (q + w X q ) 

tab b a a 

(A-51) 

The tank is rigid, therefore, the velocity of point A is sero, relative 
to the tank axis. Furthermore, the constraint condition is imposed 
when * qj,. These conditions simplify equation (4-51). Hence, 


= % (4-52) 

The impermeability constraint requires that the velocity of the liquid 
normal to the constraint surface be zero. This can be expressed as. 
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n • q, ■ 0 

D 


(4-53) 


where 

n " the outward normal vector to the constraint surface 

The normal vector is defined as the gradient of the constraint sur- 
face, 4, evaluated at the location of the liquid mass point. Thus, 


k" 
x.y.e 

(4-54) 

The tiire derivative of the constraint equation, (4-53), is required for 
the solution of the equations of motion. Therefore, 

^ (V4*qjj) ■ + (u)XV(ji) • qj^ + V4*(5 X^j^) ■ 0 

(4-55) 

From vector algebra, the following relationships are evident, 

• • 

V4* (u X qj^) - (V4 X 

(4-5ha) 


.|i 



3x 

x.y.x 3y 

x,y,z 


and 


V4Xu - -wXV4 


(4-56b) 


Thus, equation (4-55) simplifies to. 





0 


(4-57) 


This expression is equivalent to equation (4-39) and shows that the 
auxiliary equations are not functions of the velocities of the ET. 

The BT tank is axisymmetric. Therefoie, the definitions of the con- 
straint surface and fb1 matrices follow that of the two-dimensional de- 
velopment; this is also true for the baffles region of the constraint 
surface. Therefore, 
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1-0 


(4-58) 


and rcfarrlng Co Flgura 3-8 » eha bafflaa raglon 4 ls» 



0 


(4-59) 


wtier« X, y, s are mr jured in Che tank reference frame. The fb] 
matrices for Che bare and bafflea regions are. 



lb] 
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3z 
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(4-60) 


(4-61) 
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where 


"'V " 

■ Che coordinates of Che i Ch liquid mass 
point 

xj. yj, - 

■ Che velocity components of the 1 Lli liquid 
mass point 
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The tine derivative* of the (b1 matrices are, 


[b] (U 
J J 



~^3 '''' 

(Rq - Zj) ^ 



(4-62) 


and 


Ib] 


J 




(4-63) 


where 


m 

k is defined as the Zj velocity coefficient in equation (3-3S), 
Chapter 3.0. 
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4.5 THREE-DIMENSIONAL INCOMPRESSIBILITY CONSTRAINT 


The three-dinensionel liquid incompressibility wss modeled with s solid 
parallelepiped element. This element (also used to define viscous 
shear forces) maintained an average incompressibility within its vol- 
ume; the volume of any fraction of the eleaient could vary, but the ele- 
ment as a whole maintained constant volusie. This averaging technique 
was accomplished by dividing the parallelepiped into two unique ar 
rangements of tetrahedra as shorn in Figure 4-5. The voluaw of each 
arrangement (consisting of five tetrahedra) equals the volume of the 
parallelepiped. The parallelepiped volume was defined as one-half the 
sum of the ten tetrahedra volumes. The incompressibility constraint 
coefficients, [b] , for a given parallelepiped, are formed for each of 
the two unique tetrahedra arrangements. The constituents of the two 
[b] matrices are then averaged, term by term, to arrive at the final 
incompressibility constraint matrix for that parallelepiped element. 

The incompressibility constraint coefficients for each tetrahedron are 
derived by requiring the volume of the tetrahedron to remain constant 
The tetrahedron volume can be written as a scalar triple product. Re- 
ferring to Figure 4-6, the volume for a tetrahedron with the local ele- 
ment node numbering Z , m, n, o, can be written as, 

¥ = 1/6 A*(B X C) * Constant (4-64) 


The rate of change of the tetrahedron volume must equal zero. 
Therefore, d « 
dt 

+ A^CB^C^ - B^C^) + A2(B^C^ - B^C^) + A2(b3C^ - B^C^) 

+ A3(B^C2 - B^C^) + A2(Bj^C2 - B^C^) + A3(B^C2 - B 2 C 2 ) " 0 


Collecting the terms of equation (4-65) by the node point velocities of 
the tetrahedron, the components of the fbl matrix for that particular 
tetrahedron are formed. 

Node £ Components 


[(B 3 C 2 - B^Cj) - (A 3 C 2 - 
[(B2C3 - B^C^) - (Aj^C3 - 
[(B 2 C 1 - B^C2) - (A2C^ - 


AjCj) - CA2B3 - A3B3) 11^ 

A 3 C 1 ) - (A 3 B 3 - A 3 B 3 >]y^ 

A1C2 ) - (Aj^B2 - A2B^)]z^ 


(4-66a) 

(4-66b) 

(4-66c) 
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(4-67a) 


(4-67b) 

‘*2<=1 - V 2 J'. 

(4-67c) 

Node n Components 


[AjBj - *3821*^ 

(4-68a) 


(4-68b) 


(4-68c) 

Node 0 Components 


[B2C3 - B3C3li„ 

(4-69a) 

‘® 3 "l - 

(4-69b) 

1^2 - 

(4-69c) 

The components of each tetrahedron (for the two arrangements) were as- 
sembled to form two global (b| matrices for the parallelepiped ele- 
ment. This assembly used the node point numbering key indicated below 


each parallelepiped in Figure 4-5. The two global (b1 matrices were 
averaged, term by term, to form the incompressibility constraint coef- 
ficients for a given parallelepiped element. The parallelepiped [bl 
matrices were assembled to form the incompressibility constraint matrix 
for the liquid. 


The time derivative of the constraint matrix is required in the defini- 
tion of the auxiliary equation (4-39). These time derivatives were de- 
rived using an averaging technique similar to that used in deriving the 
[bl matrix. Using equations (4-66) thru (4-69), the components of the 
[b] matrix for a particular tetrahedron are, 
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2 3 

Nodes: 3, 7, 8, 6 


2 3 

Nodes; 3, 1, 6, 8 


Figure 4-6. Tetrahedron Definition of Parallelepiped Element 
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^ ■ V^" ‘ 

- Aj^l" + A2j” + A^k" 

B - (x - X )i” + (y - y.)J" + (z - 
n X, m 1C n 

- B^i*' + B2J" + B^k" 

C - (x^ - x,)i” + (y^ - y,)J” + (z - 


x,)k- 


n 


n 


z,)k- 


z,)k- 


Cj^l” + C^j" + Cjk” 


Figui'e 4-6. Tetrahedron Vector Definition 

Node Components 

l"®2^3 + - B2CJ + BjC 2 - AjC 2 + A2C2 - 

AjCz + A2CJ - A2BJ + AjB 2 - A2BJ + Ajij!*^ 


[-BjCj + B2C3 - B3C2 + BjCj - A3C3 * A3C3 - 
A3C3 t A3C3 - A3B3 + A3B3 - A3B3 + A3S3I;, 


[-B1C2 + B2C3 - B3C2 + B2C2 - A2C3 + A3C2 - 

I 

A2C3 + A3C2 - A3B2 + i2®l - 
Node m Components 

[A^C^ - A2C3 + A^C2 - A2C3lx^ 

[A1C3 - A^Cj^ + Aj^C^ - AjCj^ly^ 

[A2C1 - Aj^C 2 + A2 Cj^ - Aj^C2lz^ 

Node n Component ■ 

‘V 3 - *3“2 + *2®3 - *3”2>*„ 

‘*3“l - *l“3 * *3“l - 

1*1®2 - ‘2»1 ^ V 2 - 

Node o Components 

[B 2 C 3 - ^C 2 + B 2 C 3 - B 3 C 2 i;^ 


(4-70«) 


(4-70b) 


(4-70c> 


(4-7U) 

(4-71b) 

(4-71c) 

(4-72e) 

(4-72b) 

(4-72c) 


(4-73«) 
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FLjure 4-7. Definition of Stress Tensor for Liquid 


These components were evereged and assembled to form the incompressi- 
bility constraint matrix [b] for the liquid. 


4.6 DERIVATION OF THE LIQUID SHEAR F0IU:ES IN THREE DIMENSIONS 


A three-dimensional parallelepiped finite element was used to represent 
the shear forces within the liquid. This element was derived from the 
liquid stress-strain rate relationships using a variational mechanics 
approach. 

Consider the stresses on a differential volume of site Ax, Ay, Az (Fig- 
ure 4-7). Furthermore, assume the stress field results from liquid re- 
sistance to strain rate (deviatoric stresses) and the mean pressure on 
the differential volume. 


The stress-strain rate relationships for a viscous liquid have been de- 
rived [Ref 15]. Hence, the stresses on the differential volume can be 
expressed in terms of the strain rates as, 


0 

XX 


« -p 


a 

yy 


« -p 


9u 

(4-75a) 

3x 


3v 

(4-75b) 

3y 



4-26 




where 

u, V, w ■ the velociciee of 
tione , 
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(4-75c) 

(4-75d) 

(4-75e) 

(4-75f) 

the liquid in the tank x, y, a direc- 


a “ coefficient of viacoeity, 

P ■ mean preaaure. 


3u 3v 3w 

3x 3z 3z 


The velocity-strain rate expressions can be written using tensor nota- 
tion. 


e 


ij 


2 


(4-76) 


Furthetaore, the liquid has been assuaed incoapressible. Hence, 


3u 3v 3w 

3x 3y 3z 


0 


(4-77) 


Using equation (4-76) and the ineoapressibility statement, equation 
(4-77), the stress-strain rate relationships can be written in, tensor 
notation, 


a 


IJ 


-P 



+ 2y e 


iJ 


(4-78) 


where 


6 


iJ 


Knonecker Delta 


1 i ^ j 6-0 

I i - j 6-1 
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Defining the deviatoric atrets 

- ' - 2p E, 


Ij 


ij 


(4-79) 


the functional of the work done by the liquid atreaaea can be written 


^ " 2 ^Ij 2 ~ 

The variation of thia functional ia, 


) dVe 


(4-80) 


6J 


•-( 

4 




Ve ^Ij 2 /ve °ij “ 2 /ve 


- -f 

2 Jv 


ve "'ij 


Since only incompreaaible liquida are being conaidered, 


and 


'll* ° 


«'ii* “ 


Thua, the variation of the work integral aimplifiea to, 


6J 


2 JV 




(4-81) 

(4-82) 

(4-83) 

(4-84) 


Uaing equation (4-79), the variation of the deviatoric atreaa can be 
expreaaed aa a function of the attain rate variationa, further simpli- 
fying equation (4-84). 


5J 


■/' 




(4-85) 


This integral can be written xn matrix form as. 




6J ■ 2p|„ [(5e 6e 6t 6e 5e 6e 6c 6e 6c ] 

fVe XX xy xz yy yx yz zz zx zy 


XX 


xy 
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dVe 


zy 
(4-86) 
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The liquid etrein retee can be expreeeed in tenM of the liquid veloci- 
tiet thru the velocity-etrein rate relationshipa (equation 4-76) • Fur- 
thcmore, the velocity field within a finite eleMnt region can be ap- 
proximated by the nodal point velocitlee of the finite ei^'ient. 
wise, the atrain rate components can alao be approximated by the nc'ial 
point velocitiea. Using a tetrahedron finite element with four nodc- 
pointe, the velocity field and, her'ei the atrain rate componenta with- 
in the element, can be expreeeed ae, 


and 


“2 


(4-87) 
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3x 


3112 


3N, 

3x. 



(4-88) 


where 


N, , N 


2» ^3* *^4 


the element ehape function aasoci- 
ated with each node point; 


* 1 ’ 


^2* ^3’ ^4 • 


nodal point velocitiea. (q, 


'!• >'l- ‘l> 


Substitution of the strain rate components into the velocity-strain 
rate relationships (equation 4-76) yields, 
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(4-89) 


4-29 





ORIGINAL PAGE IS 
OF POOR QUALITY 


[311^] {q^} 
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2 3y 


(4-90) 


X, 


The variation of the strain rates can be replaced by the variation of 
the nodal point velocities. Therefore, substituting equation (4-89) 
into equation (4-86) yields. 


6J = 2p 




[3Nl’’ [3N] {q} d¥e 


(4-91) 


The nodal point viscous forces, F , for a tetrahedron can be shovm 
from the above expression to be. 



(4-92) 


where 






- T 

l3Nr 


[3N] d¥e 


4.6.1 Tetrahedron Soape Function Definition 

The specific form of the partial derivatives and the averaging of the 
dissipative forces within the parallelepiped element must be deter- 
mined. The shape function used in the partial derivatives must be at 
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least first order approximation functions. A linear approximation to 
the displacement field is, 


hj “ + xa^ + ya^ + za^ 



where 

» displacements within the tetrahedron 
- coefficients of the interpolation polynomial 


(4-93) 


A local coordinate system originating at the centroid of the tetrahe- 
dron can be developed. The nodal point displacements can then be writ- 
ten using the linear approximation function and the local nodal point 
coordinates. 



%ihere 


h 


i ■ 


X 


i’ 



displacements of the i th tetrahedron node point 
(hj^ * *1* ^i* ^1* i"l» 2, 3, 4) t 

local coordinates of the i th node point of the 
tetrahedron. 


Solving equation (4-94) for {a} , and substituting into equation (4-93) 
the element shape functions are obtained. 


h 


j 


where 



(4-95) 


hj^ “ deflection in the j th direction at node 1, etc; 

+ bj^x + c^y + d^^z 

N3 ■ + b^x + c^y + d^z 

\ ‘’4* *=4^ ^*4^ 

*These shape functions are equivalent to volume coordinates for a te- 
trahedron. Further discussion of volume coordinates can be found in 
Reference [13]. 4-31 
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The partial derivatives of the shape functions, used in equation 
(4-92), can be written as. 






1 . 2 , 3 , 4 


(4-96a) 

(4-96b) 

(4-96c) 


These partial derivatives are constants that can be placed outside the 
integral of equation (4-92). The volume of the tetrahedron is obtained 
through equation (4-64). The nodal point forces for a parallelepiped 
were determined by averaging the tetrahedron nodal point forces for the 
ten tetrahedra of each parallelepiped (Figure 4-5). These forces form 
the generalized forces Qf applied to the liquid (equation 4-38). 

4.6.2 Viscous Formulation Checkout 


The viscous force formulation was checked by calculating the forces 
that result from shearing a unit cube of liquid. Referring to Figure 
4-8a, the shear stress and resulting shear forces can be expressed as, 


du 

T = F/A = u— (4-97a) 

dy 

or 


du 

F = Am (4-97b) 

dy 

The magnitude of the force on the cube for a unit velocity and unit 
surface area equals the viscosity, p . The forces on the cube produced 
by the viscous force subroutine are shown in Figure 4-8b. The sum of 
the nodal forces in the X direction equals p, verifying the formulation 
section 4.6.1. 
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Figure 4~8a. Shear Stress on Unit Cube 



l/4p 



Figure -l-Sb. Foroc Results from Subroutine VISCOUS | 

Figui'e 4~8. Visoous Force Formulation Validation i 
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4.7 COMPUTBR PROGRAM IMPLEMBHTATION 


Numerical techniquea similar to Choaa uaad in the two-dimanaional model 
were required to tranafotm the aquation of motion into a workable com- 
puter program. The oMthoda uaad in initial liquid poaitioning and nu- 
merical integration step aiae eatimation. parallel thoae uaad in the 
two-dimensional model. The meah generation used a three-dimensional 
mapping function, which is discusaed in the subsequent paragraphs. In 
addition, specialised program logic, used to reduce computational 
costs, will be discussed. 

^•7.1 Parallelepiped Finite Element Mesh Generator (CUBIT Program) 

The liquid volume is divided into parallelepiped finite eleisents for 
the LAMPS3 isodel. The liquid volume is bounded by the tank wall and 
free surface. The tank boundary was analytically defined as a dis- 
torted ellipsoid function, 4 < The FLUDCG subroutine [Ref 1] was again 
used to obtain an approximation of the free surface geometry. The 
CUBIT program considers only axisymmetric orientations of the liquid. 
Asytmetric mesh arrangements were obtained by the techniques discussed 
in Section 3.4.2. 

Curvilinear coordinates were used to sMp cubic meshes from the C~n“V 
space into a volume in rectangular coordinate space. This voluise has 
eight faces, each having parabolic curvature. This mapping concept is 
shotm schematically in Figure 4-9. 



4~3. Three- Dimnaional Curvilinecop Happing ; 

I 


34 


Mapping is achieved through a set of "serendipity" shape funr ions 
[Ref 13} derived from a 26-noded cube. The cube has nine p Je points 
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defining each face. Thia allowa each face to be deacribed aa a quadric 
aurface. In tenta of theae ahape functiona, the napping relationahipa 
can be written aa. 


26 

X (C, n, v) - j; N. Xc 
1-1 ^ 

26 

Y (t. n, V) - r N Yc 
i-1 
26 

z (C, n , v) - >: N Zc 

1-1 


(4-98a) 

(4-98b) 

(4-98c) 


where 

- (l-v.)(l-0(l-n) - 1/2(N2+Ng+N^l - l/^fN,^+N22+N25l 

- 4 (l-v)(l-0(l-'l) - l/2{N2j+N22l 

- (l-v)(l-n )(0 - 1 / 2 [N 2 +N^+N^q 1 - l/MN 2 j^+N 22 +N 23 l 

- 4 (l-v)(l-n)(:.)('1) - 1/2IN2J+N231 

N3 - (l-v)(0(n) - 1/21N^+Ng+N^j] - 1 / 4 (N 2 j+N 23 '*-N 2^1 
Ng - 4 (l-v)d-O(O('l) - 1'21N2^+N2^1 

- (l-v)U-t,)(n) - l/2(Ng+Ng+Nj2l - 1/AIN2 j+N2^+N25] 
Ng - 4 (l-v)(l-0(l-'i)(n) - l/2(N2j+N25l 

Ng - 4 (l*“f.) (1“H ) (I- v) (\') - l/2fN22'*’^25^ 

- 4 (l-v)(l-n)(v)(0 - 1/2(N22+N23) 

- 4 (l-v)('1)(f.)(n) - 1/2[N23+N24) 

Ni 2 - 4 (l-v)d-O(vU'O - 1/2[N2^+N25l 


4 “ 



Or rCG^^ •■ 1 -klOY 

N^3 - (v)(l-0(l-n) - l/2[Ng+N^^+N2Q] 

- 1/4[N22-HN25+N2g] 

- 4 (l-C)(l-n)(v )(0 - 1 / 2 [N 22 +N 26 l 

Ni 5 - (v)(0(l-n) - l/ZtN^Q+N^^+N^^] 

- 1 / 4 [N 22 +N 23 ‘*'” 26 ^ 

* 4 (l-n)(v)(0(n) - 1 / 2 IN 23 +N 2 g] 

» (v)(0(n) - 1/2[N^^+N3g+N^gl - l/^tN 23 +N 24 +N 26 ] 

Ni8 = 4 (l- 5 )(v)(n )(0 - 1/2[N2^4-N2g] 

Nj^g = (v)(l-0(n) - 1 / 2 [N 32 +Ni 8 +N 2 o 1 

' 1 / 4 [N 2 ^+N 25 +N 2 g] 

N2q - 4 (l-n)(l-f.)(v)(n) - 1 / 2 [N 25 +N 26 ] 

N23 * 16 (l-v)(l-n)(i-0(n)(C) 

N22 - 16 (l-v)(l-0(l-'i)(0(W 
N23 * 16 (l-n)(i-v)( 4 )(n)(v) 

N2^ = 16 (l-0(l-v)(v)(n)(c) 

N25 - 16 (l-n)(l-0(l-v)(n)(v) 

N28 - 16 (l-n)(l-c)(n)(v )(0 


h- 3 (. 




O “ face center node point 


Mesh definition is obtained by defining the number of divisions along 
sides 1, 2 and 3. The ^-n-v space is then divided into the specific 
number of divisions, and the C~h~v coordinates are obtained. These co- 
ordinates are then transformed, by the shape functions, into x-y-z rec- 
tangular coordinate space. After mapping, if the mesh boundary joints 
do not match the tank surface as defined by the distorted ellipsoid 
constraint surface (due to parabolic mapping), the boundary points are 
adjusted so that they lay on the constraint surface. 


4.7.2 ”B-MATRIX" Programing Logic 

Several specialized subroutines were developed to minimize the computa- 
tional costs of the LAMPS3 program code. These routines reduced the 
number of .operations on zero-value elements in the constraint matrices, 
[b] and [b]. Subroutine logic was be ed upon storage and manipulation 
of only the terms within the bandwidth of a given row of the [b] ma- 
trix. These routines are: 


Operation Description Subroutine 

[A] ■■ [B] {V} Multiplication of BXV 

Matrix and Vector 
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[a] ■ [b]^{V} Multiplication of BTXV 

Matrix Transpose 
and Vector 



[A] ■ [C] {V} Decomposition by Gauss SKY 

Elimination and 

Forward-Backward 

Substitution 


4.8 TEST CORRELATION TO LAHPS3 MODEL 

A limited correlation analysis was performed with the three-dimensional 
model. The model was correlated to the KC-135 aircraft teat data. The 
three axes of test measured accelerations were applied to the model, 
providing three-dimensional force and moment time histories. The force 
and total impulse correlation results for KC-135 Test 1.2.2 are shown 
in Figure 4-10. The analytical Z-axis force and total impulse resemble 
the test data. A similar result is obtained for the Y-axis force and 
total impulse. The analytical X-axis force time histories also looked 
reasonable, based on X acceleration time histories. However, since 
X-axis forces were not measured, no correlation is presented for this 
axis. The analysis employed a 2x2x2 mesh (8 elements, 24 node 
points). 

The results of Figure 4-10 can be compared to those of Figure 3-26, 
the two-dimensional analytical correlations to the KC-135 test 1.2.2 
data. The two-dimensional model employed a finer mesh, containing 
32 elements. Even though the correlation In Figure 4-10 (with the 
coarse 8 element mesh) Is good, the correlation can be Improved by 
using a finer mesh. 

The three-dimensional model produces results similar to the LAMPS2 
model at reduced computational costs (due to special programming 
techniques. Section 4.7). Additional correlation analyses using 
various Integration steps and mesh sixes could be used to further 
optimize the model. However, based on the correlation performed, 
the LAMPS3 model Is a notable Improvement over the two-dimensional 
LAMPS2 formulation. 


I 
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The test data obtained in the Phase I and Phase 11 test programs have 
greatly contributed to the understanding of lo%f-g propellant reorienta- 
tion during orbiter/ET RTLS separation. Evaluation and comparison of 
the drop tower and KC-13S data has resulted in the following conclu- 
sions regarding scaling and the influence of various parameters on re- 
orientation. 

1) The external tanks' slosh and antivortex baffles have a significant 
effect on reorientation character and forces. The baffles create 
flow turbulence, which dissipates energy and reduces reorientation 
forces. 

2) Variation in tank fill volume, in general, has little effect on re- 
orientation character (in the range of IZ to 25Z volume). However, 
with baffles installed, the antivortex baffle causes small fill 
volumes (IZ to 2Z) to travel inside to the tank, instead of along 
the tank wall. 

3) In the range of Bond and Reynolds numbers expected during RTLS, and 
simulated during testing, liquid viscosity does not have a signifi- 
cant influence on the character of bulk liquid reorientation. 

4) The liquid motion observed during drop tower and KC-135 testing is 
representative of that expected during RTLS separation. Froude 
number scaling, when B0>10 and Re>30, is valid. 

The multiple mass, finite element models (LAMPS2, LAMPS3) developed 
during Phase III of the study, are a significant improvement over the 
single-point mass model currently used in the SVDS program. Correla- 
tion between the LAMPS2/LAMPS3 models and test data has resulted in the 
following observations and conclusions. 

1) The lumped mass formulation is recommended. Simulation results 
using a lumped mass matrix, correlate better to test data, and 
computational costs are significantly reduced over runs using a 
consistent mass matrix. 

2) Correlation of the models to both drop tower and KC-135 data pro- 
vides additional validation of the model scaling concept, i.e., the 
analytical models are based on the premise that Froude number scal- 
ing is valid. 

3) Force correlation is improved by reducing the integration interval 
(dt) and by using a "fine," finite element mesh. However, an ul- 
trafine mesh is not required to achieve satisfactory impulse and 
total impulse correlation. A 4x4 mesh has been found adequate, for 
the LAMPS2 model, for both drop tower and KC-135 simulations. This 
same mesh should be adequate for full-scale ET simulation. 
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4) Higher fill volumes require a finer mesh to achieve good correla- 
tion to test data. Volumes >15% will require a meah finer than the 
4x4 mesh recommended above. 

5) The three-dimensional LAMPS3 model more accurately represents the 
actual liquid mass distribution in the tank. 

6) Analytical simulations of the baffled tank do not correlate exceed- 
ingly well to baffled test data. The best correlation has been 
achieved using only a single baffle in the analytical laodel. When 
simulating the entire baffled region, the liquid gets "trapped” and 
numerical difficulties occur. Neglecting the baffles, for actual 
RTLS simulations, will result in conservative reorientation forces 
and is recommended at this time. 

Overall, the objectives of this study have been accomplished. However, 
additional work is needed to refine the analytical tools that have been 
developed. In addition, application of the models to other intact 
aborts should be considered. The following are recoianendationu for 
additional study. 

1) Testing (drop tower, KC-135) should be conducted for higher fill 
volumes to obtain data for anticipated intact aborts. 

2) Applicability of the LAMPS2/LAMPS3 models in simulating higher fill 
volume aborts should be studied. Correlation to test data should 
be used in this assessment. 

3) Additional work needs to be done in the area of baffled tank simu- 
lation. 

4) Initialization of the liquid position should be modified. The 
technique discussed in Chapters 3.0 and 4.0 has some numerical dif- 
ficulties. Mapping directly to the asymmetric initial liquid 
orientation should be investigated. 

5) Surface tensis... force effects should be incorporated into the 
models in order to allow simulation of the low Bond number regime. 


5-2 



ORIGINAL PAGE IS 
OF. POOR QUALITY 


6.0 REFERENCES 


1. Berry, Robert L. , and Tegart, James R.: Experimental Study of 

Transient Liquid Motion in Orbiting Spacecraft, Interim 
Report. NASA CR-144003, February 1975. 

2. Berry, Robert L. , and Tegart, James R. : Experimental Study of 

Transient Liquid Motion in Orbiting Spacecraft, Final Report. 
NASA CR-144213, February 1976. 

3. Berry, Robert L., and Tegart, James R.: Analysis and Test for 

Space Shuttle Propellant Dynamics (l/60th-Scale Model Test 
Results). MCR-78-523. (Contract NAS9-15302) , March 1978. 

4. Berry, Robert L. ; Tegart, James R. ; and Demchak, Leonard J.; 

Analysis and Test for Space Shuttle Propellant Dynamics 
(1/lOth-Scale Model Test Results). 2 vols. MCR-79-560. 
(Contract NAS9-15302), May 1979. 

5. Tegart, James R. : Zero-G Aircraft Test Plan. ZGTP-001 Contract 

NAS9-15302), June 1978. 

6. Williams, Robert K. : Test Procedure for External Tank Fluid Motion 

Test. NASA TTA-T-2P635 , 17 July 1978. 

7. McCarthy, J. F., Jr.: Zero-G Propulsion Problems. Jet, Rocket, 

Nuclear, Ion, and Electric Propulsion; Theory and Design. 
Edited by W.H.T. Loh. Springer-Verlag , New York, 1968. 

8. West, R. C., ed.: Handbook of Chemistry and Physics. 53rd 

ed. The Chemical Rubber Company, Cleveland, Ohio, ld72. 

9. Jasper, J. J.; The Surface Tension of Pure Liquid Compounds. 

Journal of Physical and Chemical Reference Data, Vol I, No. 4, 
1972. 

10. E. I. du Pont de Nemours & Co.; FC-114B2, Dibromotetra- 

f luoroethane . Freon Technical Bulletin, B-4B. Wilmington, 
Delaware . 

11. Coney, T. A.: Surface Tension Viscosity and Density Measurements 

of Two Fluorocarbon Solvents. NASA TMX-1862, NASA Lewis 
Research Center, Cleveland, Ohio, August I960. 

12. Stephens, D. G. : Experimental Investigations of Liquid Impact in a 

Model Propellant Tank. NASA-TND-2913 , NASA Langley Research 
Center, Hampton, Virginia, October 1965. 

13. Zienkiewicz, 0. C.: The Finite Element Method. McGraw-Hill Book 

Compduv, Maiden Head, Berkshire, England. 3rd ed. Copyright 
1977. 


6-1 



b . Bowman, T. B>: Response of the Free Surface of a Cylfndrlcally 

Contained Liquid to Off-Axis Accelerations. Proceedings of the 
1966 Heat Transfer end Fluid Mechanics Institute, Stanford 
University Press, 1966. 

15. Lamb, Horace, Sir: Hydrodynamics. 6th ed. Dover Publications. 

New York, 1932. 



APPENDIX A— VIRTUAL POWER CONCEPT 
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The virtue! pomr concept ie e convenient way to apply the principle of 
virtual work to mechanica problem. Thia appendix diacuaaes the con- 
cepta of the virtual power aMthod. Thia documntation ia an excerpt 
from the "Skylab Docking Maneuver Simulation Nathemtical Model Report 
(Final),” ED- 2002 1395, 1971, authored by Carl Bodley and Colton Park. 

The concept of a virtual diaplacement and the principle of virtual work 
provide very effective toola in the analyaia of a mechanical ayatem. 

The utility of the principle of virtual work is minly due to being 
able to eliminate the action of the inner forces necessary to mintain 
prescribed kinemtical conditions. For example, a deformable structure 
my be analyxed for stress and strain by prescribing certain poaaible 
modes of deformation. These modes are a finite nu^er of kinemtical 
conditions with which the virtual diaplacemnts of the ayatem must com- 
ply. Because the virtual diaplacemnts are in harmony with the given 
kinemtical conditions, the action of the inner forces my be elimi- 
nated. The numoer of equations, equal to the number of modes of defor- 
mtion (or degrees of freedom assumd), is amller than the number of 
particles contained in the system. 

The principle of virtual work ir a variational principle that allows a 
given mechanical system to be analysed without considering the infinity 
of particles and the infinity of inner forces acting between them. The 
tr,. 'sformatlon that relates the rectangular coordinates of the system's 
particles (X^^, Y^, Z() to the generalised coordinates (q]^, q2, qn) 
Is used to transform z\e virtual work expression that Is In terms of 
the 6X^, 6Yj^, and the resultant particle forces; 


That is 


6W 




is transformed to 


6W - ^ Qj 

j-1 


(A-1) 


(A-2) 


A-1 
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The analytical form of the generalised forces, Qj, is obtained 
through this transformation. The transfonsation that relates the par- 
tides' rectangular coordinates to the generalised coordinates is in 
essence, the prescribed kinematical conditions of constraint: 

(**1’ ^2 '’n' 


The virtual displacements j r^ may be expressed as 


6r. 

X 



6r. 


(A-4) 


thus, from (A-1) 


N 


i=l 


F 

1 




(A-5) 


and it follows that 




(A-6) 


Equation (A-6) shows the transformation relating the particle forces to 
the generalized forces. By the principle of virtr<il work, only those 
Fi that are impressed forces contribute to 6W, and therefore, to 
0;^. If the system is a flexible solid, some of the Ff may be con- 
sidered to be resultant stresses acting on virtual strains that are 
consistent with (equation A-A); the system has the capacity for storing 


A-2 





tential energy (0). Qj^ aey be considered to be composed of contribu- 
tions from ell the conservative impressed forces and from the position- 
dependent (or conservative) forces, 



15;; - Z (v“)i • 

i-1 


(A-7) 


where L is the total of the unconservative impressed forces. In de- 
fining the explicit form of the Qj^ it is necessary to form the deriv- 
atives (3rj^)/(3qj^) . Given the explicit form of the kinematical 
conditions, such as (equation A- 3), one needs only to differentiate to 
express the required derivatives. However, often in practice the form 
of the transformation (equation A-3) is only implied and the 
derivatives are difficult to obtain. It is generally possible to 
easily express the particle velocities and the time (t), simply by 
writing well known kinematical expressions that comply with (equation 
A-3); 

J) 

^2’ ^1’ ^2* **'♦ (A-8) 

also 


‘ .!!i. 

3q. 3t 

j=l ^ 


(A-9) 


The derivatives being sought appear in (equation A-9); tne 
accomplishment has been to express rj in a simple "kinematical 
fashion" and equate the result to that given by (equation A-9) 
obtaining explicit expressions for (6r^)/(6o.) v;ithout needinp 
to differentiate (equation A-3) or witnout e^en needing to know the 
explicit form of A-3. Since the second term on the right of (equation 
A-9) does not enter into the comparison, it only needs to be recognized. 

Virtual velocity can now be defined. The virtual velocity concept is 
analogous to that of virtual displacements, in that the qj^ can be 
imagined to assume arbitrary values independent of time. Particle ve- 
locities take on values consistent with the conditions 


n 


r* = 


k=l 




(A-10) 

A-3 
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The difference between (equation A-9 and equation A-10) is the term 
Orj^)/Ot), which vanishes since "time is held constant" '.*'ile 
impressing arbitrary values of qjj^ similar to expressing virtual 
displacements. Thus, tne virtual velocity of the i^^ peirticle is 
defined to be that which satisfies (equation A-10) for arbitrary values 
of the independent qjjj. The superscript (*) is used to 
distinguish between virtual velocities and real velocities. 

may be expressed from simple well known kinematical relation- 
ships and these results compared with (equation A-10) to obtain the 
explicit form of the derivatives (3r^)/(3q ), again without needing 
to know the explicit form of (equation A- 3;. 

In conjunction with virtual velocity, virtual power, 6w, is conceived 
to be analogous to virtual work. It is an expression involving virtual 
velocities rather than virtual displacements: 


N 

6 W = F. • r*. 

1 1 


i=l 

n 

/ ^ 



II 

M 


3 r .1 

(A- 11 ) 

k=l ' 
n 

II 

Xi* I-* 

7 

(A- 12 ) 

k=l 





The virtual power has served the exact same purpose that virtual work 
does, namely to provide the explicit form of the generalized forces. 
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APPENDIX B— USER GUIDE FOR LAMPS2 PROGRAM 


This appendix provides user information for the two-dimensional 
model (LAMPS2) and its mesh generator (MESH2). 

The MESH2 program is a stand-alone program, which generates the 
finite element mesh geometry based on user input data. The mesh 
data are written on a disk file, which LAMPS2 accesses directly. 
Figure B-1 shows the use and interaction of the two programs. 

'' The tank constraint surface, <ti, is developed (automatically) 

based on input tank geometry. Any axisymmetrlc tank geometry 
can be simulated by the LAMPS2 program, from a cylinder to a 
sphere. Figure B-2 shows example geometry information input to 
the programs (MESH2, LAMPS2) ^ 

t Figure B-3 shows the overlay structure of the LAMPS2 program 

I and the computer core requirements for a CDC CYBER computer. 

Table B-1 defines the input parameters for program MESH2. 

Table B-2 shows sample input data for MERH2. Table B-3 defines 
the input parameters for program LAMPS2. Table B-4 shows sam- 
ple input data for LAMPS2. 


USER INPUT 


USER INPUT 



FORCE TIME HISTORIES MOTION TIME HISTORIES 


Figure B-1. Tim- D imensional Model Program Interaction 
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IDOHET = 
OGIVE 




IDOMET = 1 
SPHERE 
TDT=RT 



RB 


IDOMEB = 1 
ELLIPSE 
RB=TR 



T 

XL 

J. 


TR 



IDOflEB = 1 
SPHERE 
TDB=RB 


IDOMET 

RT 

TDT 

lOOMEB 

RB 

TDB 

XL 

TR 


1, top dome is ellipsoidal or spherical 

2, top dome is ogive 

radius of ellipse or ogive 
height of upper dome 

1, bottom dome is ellipsoidal or spherical 

2, botton dome is ogive 

radius of ellipse or ogive 
height of lower dome 
length of tank cylindrical section 
tank radius 


Figure F-2. Tank Constraint Surface Geometry Definition 
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LAMPS2 MESH GENERATOR PROGRAM 
PROGRAM FORMS FINITE ELEMENT MESH FOR LAMPS2 
DEVELOPED BY : R.L. BERRY AND L.d DEMCHAK 

FOR NAS9- 15302 

LAST REVISION APRIL 1980 

TANK VARIABLES EXPECTED. 


FOR THE TOP DOME: 

IDOMET = 1, TOP ELLIPSOIDAL OR SPHERICAL DOME 
- 2 , TOP OGIVE DOME 
RT = RADIUS OF ELLIPSE OR OGIVE 

TDT = HEIGHT OF ELLIPSE OR OGIVE (TRUNCATED HGT OF OGIVE) 

FOR THE BOTTOM DOME: 

IDOMEB = 1, BOTTOM ELLIPSOIDAL OR SPHERICAL DOME 
= 2. BOTTOM OGIVE DOME 
RB = RADIUS OF ELLIPSE OR OGIVE 

TDB = HEIGHT OF ELLIPSE OR OGIVE (TRUNCATED HGT OF OGIVE) 

XL = LENGTH OF PROPELLANT TANK CYLINDRICAL SECTION. 

(L UNITS) 

TR = TANK RADIUS. (L UNITS) 

FLUID POSITION AND GRID VARIABLES EXPECTED 


THETAX = ORIENTATION OF TANK (DEG), IF .NE. 0 ITERATE TO 
POSITION DEFINED 

NSIDE1 s NUMBER OF DIVISIONS OF SIDE 1 OF FLUID GRID, 

IF NSIDE1 .LE. 0 READ IN PREVIOUSLY FORMED GRID DATA 
NSIDE2 = NUMBER OF DIVISIONS OF SIDE 2 OF FLUID GRID 
NMESH = 1, SIDE1 ,SIDE2, SIDES ON BOUNDARY, SIDE4 ON 
FREE SURFACE 

NMESH = 2, SIDE1,SIDE4 ON BOUNDARY, SIDE2.SIDE3 ON 
FREE SURFACE 

MESH GENERATION CAPABILITY 


THE PROGRAM IS CURRENTLY SET UP FOR A MAXIMUM 
GENERATION OF 81 NODE POINTS. THIS IS EQUATED 
AS FOLLOWS: NUMNP = (NSIDE1 + 1) X (NSIDE2 + 1) 

»*NOTE*» THE LAMPS2 PROGRAM HAS ITS OWN CRITERIA 
FOR THE MAXIMUM NUMBER OF NODE POINTS. THIS IS 
CURRENTLY SET AT NUMNP = 49 FOR COST SAVINGS. 

FLUID PARAMETERS EXPECTED 


FDEN = MASS DENSITY OF FLUID 
PCVOL = PERCENT TANK FILL. LE. 100. 

AMTOT = ACTUAL LIQUID MASS 

PLOT TITLE INFORMATION EXPECTED 

PTIT = ALPHA-NUMERIC PLOT TITLE (4A10) 

Table B-1. Program MESH2 Input Parameter Definition 
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BAFFZLs 

BAFFZUs 

BAFHGT= 

BAFCON= 
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BAFFLE(S) PARAMETERS EXPECTED 


O.NO BAFFLES ADDITIONAL INPUT NOT REQUIRED 

1.0NE BAFFLE SET BAFFZL=BAFFZU, INPUT BAFCON 

N,N NUMBER OF BAFFLES , BAFCON OVERRIDDEN 

LOWER BAFFLE LOCATION 

UPPER BAFFLE LOCATION 

EFFECTIVE HEIGHT OF BAFFLE 

EFFECTIVE WIDTH OF BAFFLE 


able B-1 


(aOTial) 
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MESHGN S.W.WIRTH 

MESH GENERATION FOR CASE 19 5 % FILL;MESH TYPE 1, 5X7 

ONE BAFFLE IS UTILIZED FOR THE FULL BAFFLE DISTANCE 
MESH/MASS PRG 5X7 MESH VERSION 
PSINPUT 
XL=1.31, 

TR=2.75, 

RT=10. 16, 

TDT=6.46, 

ID0MET=2, 

RB=2.75, 

TDB=2.06, 

THETAX=0.0, 

PCVOL=5. , 

ID0MEB=1, 

FDEN=2.0213E-4, 

AMT0T=1 .5666E-03, 

NSIDE1=5, 

NSIDE2=7, 

NMESHsI . 

NBAFFsl , 

BAFFZL=1 .00, 

BAFFZU=1 .00, 

BAFHGT=0.5, 

BAFC0N=3.5, 

$ 

STOP 


Table B-2. Sample Input Data for Program MESH2 
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4M«I m «*»■«> m * I«c *** 4i 4I * If * *4i m 4< 4> * 4>4I * 4> * 4> * 4t «**«** *1^ *»4I* * 41 4> 4I «**■•<* iK XI *** 4t 

LAMPS2 PROGRAM 

> LARGE AMLPITUDE 2 DIMENSIONAL SLOSH ANALYSIS PROGRAM 

■I DEVELOPED BY : R.L. BERRY AND L.d. DEMCHAK 

1 FOR NAS9- 15302 

f LAST REVISION MAY 1980 

41 * Xc 4i « « Xt Xi X> XI * * « * XI « 4i * * * * Xi X> « * * * * Xi XI « * * * * 4> Xi * * * * * 4t XI * * « * * X> * * Xi XI Xi 41 « « * 4i » 

TAPE1 RESERVED FOR FORCE TIME HISTORY STORAGE 
TAPE2 RESERVED FOR INPUT OF EXTERNALLY GENERATED MESH DATA OR 
RESTART 

TAPES RESERVED FOR LIQUID MOTION TIME HISTORY STORAGE 
; TAPE4 RESERVED FOR RESTART DATA OUTPUT 

TANK VARIABLES EXPECTED. 


IDOMET = 1, TOP ELLIPSOIDAL OR SPHERICAL DOME 
= 2. TOP OGIVE DOME 
RT = RADIUS OF ELLIPSE OR JGIVE 

TDT = HEIGHT OF ELLIPSE OR OGIVE (TRUNCATED HGT OF OGIVE) 

IDOMEB = 1, BOTTOM ELLIPSOIDAL OR SPHERICAL DOME 
IDOMEB = 2, BOTTOM OGIVE DOME 
I RB = RADIUS OF ELLIPSE OR OGIVE 

TDB = HEIGHT OF ELLIPSE OR OGIVE (TRUNCATED HGT OF OGIVE) 

XL = LENGTH OF PROPELLANT TANK CYLINDRICAL SECTION. 

(L UNITS) 

TR = TANK RADIUS. (L UNITS) 

FLUID POSITION AND GRID VARIABLES EXPECTED 


THETAX = ORIENTATION OF TANK (DEG), IF .NE. 0 ITERATE TO 
POSITION DEFINED 

NTHETA = NUMBER OF INCREMENTS OF THETA TO ACHIEVE 
CONVERGENCE 

CONVKE = CHANGE IN KINETIC ENERGY FOR POSITION CONVERGENCE 
CONVDT = INTERGRATION STEP SIZE FOR INITIAL POSITION 
CONVERGENCE 

G1 s ACCELERATION MAGNITUDE FOR INITIAL POSITION 
CONVERGENCE 

NITERP = NUMBER OF ITERATIONS TO ESTABLISH FLUID INITIAL 
POSITON, THETAX .NE. 0. FOR ITH THETA INCREMENT 

FLUID PARAMETERS EXPECTED 


FDEN = MASS DENSITY OF FLUID 

XMUC = FLUID VISCOSITY (CENTIPOISE) , IF .LE. 0 NO VISCOUS 
FORCES 

FGLUE = GLUE FACTOR FOR FLUID TO REMAIN ON WALL, 

RECOMMEND 1.0 

Table B~3. Program LAMPS 2 Input Parameter Definition 
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PCVOL = PERCENT TANK FILL. LE. 100. 

AMTOT = ACTUAL LIQUID MASS 

INTERGRATION AND PRINT/PLOT FLAGS EXPECTED 


DT - MAXIMUM DELTA TIME INCREMENT 
TF = END TIME OF SIMULATION 
NPRNT = PRINT COUNTER, PRINT EVERY NPRNT TIME 
NPLOT = PLOT COUNTER, PLOT EVERY NPLOT TIME 
IFSAVE * SAVE STATE ,ETC. TO JUMP START LATER (TAPE4) 
IFJUMP = RESTART PROGRAM FORM PREVIOUS RUN 

(SEE SUBROUTINE dUMPRD FOR INPUT REQUIREMENTS) 

PLOT TITLE INFORMATION EXPECTED 


PTIT = ALPHA-NUMERIC PLOT TITLE (4A10) 
BAFFLE(S) PARAMETERS EXPECTED 


NBAFFs 0,NO BAFFLES ADDITIONAL INPUT NOT REQUIRED 

= 1,ONE BAFFLE, SET BAFFZL*BAFFZU, INPUT BAFCON 
- N,N NUMBER OF BAFFLES, BAFCON OVERRIDDEN 
BAFFZL= LOWER BAFFLE LOCATION 
BAFFZU= UPPER BAFFLE LOCATION 
BAFHGTs EFFECTIVE HEIGHT OF BAFFLE 
BAFCON= EFFECTIVE WIDTH OF BAFFLE 

ACCELERATION DATA EXPECTED 


ACCEL = TABLE OF TANK ACCELERATIONS: 

COL 1- TIME, COL 2- AY, COL 3- AZ 

MESH DATA EXPECTED 

TABLES DEFINING FINITE ELEMENT MESH FROM PROGRAM MESH2 ON FILE 2 
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LAMPS2 S.W.WIRTH 
TEST 19 CASE 1 LAMPS RUN 
5 % FILL WITH ONE BAFFLE 
TEST 19 CASE1--5% FILL LAMPS RUN. 
P$INPUT 
XL=1.31, 

TRs2.75, 

RT=10.16. 

TDT=6.46, 

IDOMET=2, 

RB*2.75, 

TDB=2.06, 

THETAX=13.0, 

PCVOL=5. , 

IDOMEB^I, 

FDEN=2.0213E-04, 

XMUC=1 .0, 

SMASSs4.751E-03, 

AMT0T=1.5666E-03, 

ZCGS’=-35.9045, 

YCGSsl2.7403, 

DT=0.02, 

TFsI.OO, 

FGLUEsI .0, 

NPRNT=4, 

NPL0T=4, 

NITERP*40, 

NBAFFsl , 

BAFFZLsI.O, 

BAFFZUsl.O, 

BAFHGT=0.5, 

BAFCON=3.5, 


G1=190.0, 

CONVKE=5.0E-04, 

CONVDT=0.025, 

$ 




TIME 


AY 


AZ 


ACCEL 

3 

3 



1 

1 

0.0 

12.7403 

386.07 

2 

1 

0.07 

12.7403 

•35.1320 

3 

4 

( 

1 .6 

12.7403 

-35. 1320 

0000000000 




YZCORD 

0 

2MESHGN 1 



NTYPE 

0 

2MESHGN 1 



NP 

0 

2MESHGN f 



NLINE 

0 

2MESK3N } 

DATA FROM PROGRAM MESH2 


AMASS 

0 

2MESHGN \ 



AREA 

0 

2MESHGN 



HQT 

0 

2MESHGN ' 



STOP 






Table B’-4. 


Sample Input Data for Program LAMPS2 
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APPENDIX C— USER GUIDE FOR LAMPS3 PROGRAM 

This appendix provides user information for the three-dimensional 
model (LAMPS3) and its mesh generator (CUBIT). 

The CUBIT program is a stand-alone program, which generates the 
finite element mesh geometry based on user input data. The mesh 
data are written on a disk file, which LAMPS3 acri<sses directly. 
Figure C-1 shows the use and interaction of the two programs. 

The tank constraint surface, (^, is developed ( .^tomatically) 
based on input tank geometry. Any axlsymraetrlc tank geometry 
can be simulated by the LAMPS3 program, from a cylinder to a 
sphere. The tank geometry information required is identical 
to that for the two-dimensional model (see Figure B-'^). 

Figure C-2 showu the overlay structure the LAMPS3 '-rogram and 
the computer core requirements for a CDC CYBER computer. Table 
C-1 defines the input p rameters for program CUBIT. Table C-2 
shows sample input data for CUBIT. Table C-3 defines the input 
parameters for program LAMPS3. Table C-4 shows sample input 
data for LAMPS3. 



Figure C-1. Three- Dimnaional Model Program Interaction 
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Figure C~2. LAMPS 3 
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CUBIT MESH GENERATOR PROGRAM 
PROGRAM FORMS FINITE ELEMENT MESH FOR LAMPS3 
DEVELOPED BY : R.L. BERRY, L.J DEMCHAK, A. COVINGTON 
FOR NAS9- 15302 
LAST REVISION SEPT 1980 

TANK VARIABLES EXPECTED. 


IDOMET = 1, TOP ELLIPSOIDAL OR SPHERICAL DOME 
=2, TOP OGIVE DOME 
RT = RADIUS OF ELLIPSE OR OGIVE 

TDT = HEIGHT OF ELLIPSE OR OGIVE (TRUNCATED HGT OF OGIVE) 

IDOMEB = 1. BOTTOM ELLIPSOIDAL OR SPHERICAL DOME 
IDOMEB = 2, BOTTOM OGIVE DOME 

RB = RADIUS OF ELLIPSE OR OGIVE 

TDB = HEIGHT OF ELLIPSE OR OGIVE (TRUNCATED HGT OF OGIVE) 

XL = LENGTH OF PROPELLANT TANK CYLINDRICAL SECTION. 

(L UNITS) 

TR = TANK RADIUS. (L UNITS) 

GRID VARIABLES EXPECTED 


NSIDE1 = NUMBER OF DIVISIONS ON X AXIS (TOP SURFACE) 

NSIDE2 = NUMBER OF DIVISIONS ON Y AXIS (TOP SURFACE) 

NSIDE3 s NUMBER OF DIVISIONS ON Z AXIS (DEPTH OF TANK) 

FLUID PARAMETERS EXPECTED 


FDEN = MASS DENSITY OF FLUID 
PCVOL = PERCENT TANK FILL. LE. 100. 
AMTOT = ACTUAL LIQUID MASS 

PLOT INFORMATION EXPECTED 


PTIT = ALPHA-NUMERIC PLOT TITLE (8A10) 

IFd = 1, PLOT NODE NUMBERS 
NVW = 1, PLOT FULL TETRAHEDRA ELEMENTS 
= 2. PLOT ONLY CUBES 
IFPT s PRODUCE DIAGNOSTIC PRINT OUTPUT 
IFCORN= INPUT OPTION FOR GENERAL MESH GENERATION, 

REQUIRES ADDITIONAL INPUT OF "CORNER" NODES BY 
CALL READ (3 TIMES) 


Table C-1. Prog.''an CUBIT Input Parameter Definition 
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ORIGINAL PAQf It 
OF POOR QUALITY 


BAFFLE! S) PARAMETERS EXPECTED 


NBAFFs O.NO BAFFLES ADDITIONAL INPUT NOT REQUIRED 
= 1.0NE BAFFLE SET BAFFZL=BAFFZU , INPUT BAFCON 
s N.N NUMBER OF RAFFLES , BAFCON OVERRIDDEN 
BAFFZLs LOteER BAFFLE LOCATION 
BAFFZUs UPPPR BAFFLE LOCATION 
BAFHGTs EFFECTIVE HEIGHT OF BAFFLE 
BAFCONs EFFECTIVE WIDTH OF BAFFLE 


C-4 


Table C-1. 


(ooncl) 
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MESH3D S.W.WIRTH 

MESH GENERATION FOR 3D MODEL -10 % FILL; 1 BAFFLE 
2X2X2 MESH TYPE FOR P4LMP3B 
3DMESH FOR 10% FILL; 2X2X2; 1 BAFFLE 
PSINPUT 
XL=1.31, 

TR=2.75, 

RT=10.16, 

TDT=6.46, 

IDOMET=2, 

RB=2.75, 

TDB=2.06, 

PCVOL=10. , 

IDOMEB=1, 

FDEN=2.0213E-4. 

AMTOT=3. 1212E-03, 

NSIDE1=2, 

NSIDE2=2. 

NSIDE3=2, 

NBAFF=1, 

IFCORN=0, 

IFPT=0, 

NVW=2, 

BAFFZL=1 .00, 

BAFFZU=1 .00, 

BAFHGT=0.5, 

BAFCON=3.5, 

$ 

STOP 


Table C-2. Sample Input Data for Program CUBIT 
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LAMPS3 PROGRAM 

LARGE AMLPITUDE 3 DIMENSIONAL SLOSH ANALYSIS PROGRAM 
DEVELOPED BY : R.L. BERRY AND L.J. DEMCHAK 
FOR NAS9- 15302 
LAST REVISIOr DEC 20 1980 

4c It: * * 4: * * « 4c 41 4c ♦ 4 e 4: # 4c 4: 4: 4c 4c lie * * « * 4 e * ♦ * 4( * * 4< ♦ ♦ 4c 4c * * * * « « itc » 4i 4c 4c * * * 4c * >*< * * 4= 4^ * 

TAPE1 RESERVED FOR FORCE TIME HISTORY STORAGE 

TAPE2 RESERVED FOR INPUT OF EXTERNALLY GENERATED MESH DATA 

OR RESTART 

TAPE3 RESERVED FOR LIQUID MOTION TIME HISTORY STORAGE 
TAPE4 RESERVED FOR RESTART DATA OUTPUT 

TANK VARIABLES EXPECTED. 


IDOMET = 1, TOP ELLIPSOIDAL OR SPHERICAL DOME 
=2, TOP OGIVE DOME 
RT = RADIUS OF ELLIPSE OR OGIVE 

TDT = HEIGHT OF ELLIPSE OR OGIVE (TRUNCATED HGT OF OGIVE) 

IDOMEB = 1, BOTTOM ELLIPSOIDAL OR SPHERICAL DOME 
IDOMEB = 2 , BOTTOM OGIVE DOME 

RB = RADIUS OF ELLIPSE OR OGIVE 

TDB = HEIGHT OF ELLIPSE OR OGIVE (TRUNCATED HGT OF OGIVE) 

XL = LENGTH OF PROPELLANT TANK CYLINDRICAL SECTION. 

(L UNITS) 

TR = TANK RADIUS. (L UNITS) 

FLUID POSITION AND GRID VARIABLES EXPECTED 


THETAX = ORIENTATION OF TANK (DEG), IF .NE. 0 ITERATE TO 
POSITION 

NTHETA = NUMBER OF INCREMENTS OF THETA TO ACHIEVE 
CONVERGENCE 

CGNVKE = change IN KINETIC ENERGY FOR POSITION 
CONVERGENCE 

CONVDT = INTERGRATION STEP SIZE FOR INITIAL POSITION 
CONVERGENCE 

G1 s ACCELERATION MAGNITUDE FOR INITIAL POSITION 
CONVERGENCE 

NITERP = NUMBER OF ITERATIONS TO ESTABLISH FLUID INITIAL 
POSITION, THETAX .NE. 0 FOR ITH THETA INCREMENT 
THETAX .NE. 0. FOR ITH THETA INCREMENT 

FLUID PARAMETERS EXPECTED 


XMUC = FLUID VISCOSITY (CENTIPOISE) , IF .LE. 0 
NO VISCOUS FORCES 

Table C-3. Program LAMFS3 Input Parameter Definition 
C-6 
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FGLUE = GLUE FACTOR FOR FLUID TO REMAIN ON WALL 
RECOMMEND 1.0 

AMTOT = ACTUAL LIQUID MASS 

INTERGRATION AND PRINT/PLOT FLAGS EXPECTED 


DTMAX = MAXIMUM DELTA TIME INCREMENT 
TF = END TIME OF SIMULATION 
NPRT = PRINT COUNTER. PRINT EVERY NPRT TIME 
IFPLT = DIAGONSTI FOR PLOT ROUTINES 
IFSAVE = SAVE STATE .ETC. TO JUMP START LATER (TAPE4) 
IFJUMP = RESTART PROGRAM FORM PREVIOUS RUN 

(SEE SUBROUTINE JUMPRD FOR INPUT REQUIREMENTS) 
NCHECK = 1. MINIMUM OUTPUT (DEFAULT VALUE) 

= 2, ST TE OUTPUT/ B BD ETC. . . 

= 3. ALL OUTPUT, BRING WHEEL BARREL 

PLOT TITLE INFORMATION EXPECTED 


PTIT = ALPHA-NUMERIC PLOT TITLE (4A10) 
NPLOT = PLOT COUNTER, PLOT EVERY NPLOT TIME 
NSIDE1 = NUMBER OF DIVISION ALONG X AXIS 
NSIDE2 = NUMBER OF DIVISION ALONG Y AXIS 
NSIDE3 s NUMBER OF DIVISION ALONG Z AXIS 
IPLOPT = 0, FREE SURFACE PLOT OPTION 
IPLOPT = 1, BRICK PLOT OPTION 

BAFFLE(S) PARAMETERS EXPECTED 


NBAFFs 0,NO BAFFLES ADDITIONAL INPUT NOT REQUIRED 

= i.ONE BAFFLE, SET BAFFZl.sBAFFZU , INPUT BAFCON 
= N,N NUMBER OF BAFFLES, BAFCON OVERRIDDEN 
BAFFZLs LOWER BAFFLE LOCATION 
BAFFZUs UPPER BAFFLE LOCATION 
BAFHGTs EFFECTIVE HEIGHT OF BAFFLE 
BAFCON= EFFECTIVE WIDTH OF BAFFLE 

ACCELERATION DATA EXPECTED 


ACCEL = TABLE OF TANK ACCELERATIONS: 

COL 1- TIME, COL 2- AX, COL 3- AY, COL 4- AZ 

MESH DATA EXPECTED 


TABLES DEFINING FINITE ELEMENT MESH FROM PR0GR.AM CUBIT ON FILE 2 


Table C-3, (aonol) 
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LAMPS3 S.W.WIRTH 
P4LMP3B RUN FOR 10X FILL 1 BAFFLE 
2X2X2 MESH--VISCOUS SUBROUTINE INCLUDED. DT= . 005 
LAMPS 3D FOR 10% FILL; 1 BAFF;DTs.005 
PSINPUT 
XLsI .31 . 

TRs2.75, 

RTslO.16, 

TDT=6.46, 

ID0MET=2, 

RB=2.75, 

TDB=2.06, 

THETAX=13.0, 

IPLOPTsI, 

IDOMEB=1, 

IFPLT=0, 

XMUC=0.75, 

SMASS=5.77387E-03, 

AMTOT=3. 1212E-03, 

Z0FF=35.9425. 

NTHETAsI, 

DTMAX=0.005, 

TF=1 . 10, 

PQLUE=1 .0, 

NPRT=4, 

NPL0T=4, 

NITERP=40, 

IFAPP=1 , 

NBAFF = 1 , 

BAFFZL=1 .0, 

BAFFZU=1 .0, 

BAFHGT=0.5. 

BAFC0N=3.5, 

G1=190.0, 

CONVKEsI .OE-04, 

CONVDT=0.005, 

IFSAVE=0, 

IFdUMP=0, 

NSIDE1=2, 

NSIDE2=2, 

NSIDE3=2, 


Table C-4. Sample Input Data for Program LAMPS3 
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A TIME 

^AX 

AY ^ AZ 

\ 

$ 

ACCEL 

3 

4 f 


▼ 



1 

1 

0.0 


0.0 

0.0 

386.07 

2 

1 

0.07 


0.0 

11.582 

-34.94 

3 

1 

2.0 


0.0 

11.582 

-34.94 

0000000000 






XYZCOR 

0 

-2MESH3D 





NTYPE 

0 

-2MESH3D 





NP 

0 

-2MESH3D 

! 




MEND 

0 

-2MESH3D 


DATA FROM PROGRAM CUBIT 


IDIAG 

0 

-2MESH3D 

1 




AMASS 

0 

-2MESH3D 





ELVOL 

0 

-2MESH3D 

f 





STOP 


Table C-4. (aonol) 
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